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iSTRACT  (Continuo  on  rovotmo  oldo  If  nofhprr  IdoniHy  hr  block  nirmhot) 


'^The  ability  to  predict  and  simulate  acoeleratton-time  histories  at  the  actual  mounting  locations 
of  weapon  system  components  will  provide  vital  knowledge  for  protecting  and  hardening  these 
components  to  meet  present  and  changing  mission  requirements  of  naval  ships  and  submarinea. 
Mechanical  Impedance/mobility  techniques,  long  a tool  of  the  experimenter,  have  now  evolved 
to  practical  methods  under  engineering  conditions.  This  report  summarizes  cunent  activities  in 

the  measurement  of  structural  dynamic  functions,  usable  methods  for  predicting  respoiue  to 
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20.  Abstract  (Continued) 

environmental  threats,  and  a new  method  for  environment-simulation  testing. 

llie  report  is  replete  with  example  applications  that  demonstrate  the  practicality  of  impedance/ 
mobility  techniques.  Problems  that  may  become  more  tractable  include  the  scaling  of  previous 
test  results  from  one  ship  to  another,  and  from  one  kind  of  equipment  to  another.  The  techniques 
could  be  applied  to  quasi-linear  systems,  nonlinear  systems,  parametric  identification  multivariate 
regression  analysis,  and  pulsed  environmental  simulation.  Some  projects  that  could  yield  important 
information  are  suggested  for  review  extension,  and  elaboration  by  the  naval  community. 
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PREFACE 


This  report  presents  the  results  of  a study  on  impedance-based 
motion  prediction,  scaling,  and  environmental  simulation  for  shock  applica- 
tions. G.J.  O'Hara  was  project  manager  for  the  Naval  Research  Laboratory 
(NRL);  F.B.  Safford  was  project  engineer  for  Agbabian  Associates  (AA)  and 
author  of  this  report. 

Dr.  R.O.  Belsheim,  NRL,  (now  with  NKF  Engineering)  furnished 
information  about  the  effects  of  equipment  loading  on  structural  motions. 

Or.  Vernon  Neubert,  Pennsylvania  State  University,  provided  considerable 
background  information  on  the  U.S.  Navy  impedance  projects,  as  did  R.E.  Walker, 
Waterways  Experiment  Station,  on  the  U.S.  Army  Corps  of  Engineers  projects. 

Dr.  Robert  Plunkett,  University  of  Minnesota,  served  as  a consul- 
tant during  the  development  and  evaluation  of  techniques  for  response  pre- 
diction and  pulse  simulation.  Optimization  investigations  were  made  in 
consultation  with  Dr.  George  A.  Bekey  and  Dr.  S.F.  Masri,  University  of 
Southern  California.  Transient  air-blast  scaling  analysis  was  furnished  by 
Dr.  Arthur  R.  Maddox,  Naval  Weapons  Center,  China  Lake,  California. 
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SECTION  1 


INTRODUCTION 


I. I SCOPE 

This  report  sunmarizes  the  state  of  the  art  (with  emphasis  on  recent 
advances)  for  impedance,  response  predictions  therefrom,  and  environmental 
simulation  by  pulse  techniques.  Useful  applications  of  these  methods  to 
ongoing  naval  needs  are  also  presented,  together  with  recommendations  for 
future  work  in  these  areas. 

The  background  of  impedance,  based  on  the  contributions  of  the 
U.S.  Navy,  aerospace  and  machine  tool  industries,  the  Defense  Nuclear  Agency, 
and  the  Corps  of  Engineers,  is  reviewed.  Although  impedance  is  based  on  linear 
theory,  various  feasible  methods  have  been  considered  for  projecting  the 
technique  to  parameter  identification  for  nonlinear  structural  conditions 
and  for  approximate  methods  such  as  multivariate  regression  analysis. 

Excellent  response  predictions  in  both  the  frequency  and  time  domains 
can  be  accomplished  with  impedance  records  of  magnitude  and  phase,  and  their 
associated  impulse  functions.  The  ability  to  make  time-domain  predictions  for 
large  physical  systems  is  a relatively  recent  accomplishment,  obtained  from  con- 
tinuous improvement  in  measurement  systems  and  in  digital  processing  methods. 
Response  motions  can  be  conveniently  predicted  from  impedance  and  transfer 
functions  over  a bandwidth  of  a few  hertz  to  5000  hertz.  Predictions  that  were 
compared  against  actual  records  from  field  tests  were  in  good  agreement. 

Some  informal,  but  incomplete,  comparisons  have  been  made  between 
Impedance-based  predictions  and  finite  element  models.  These  comparisons  were 
not  satisfactory  because  of  differences  in  test  ranges  and  finite  element  model 
assumptions.  A more  formal  comparative  study  is  now  being  made  of  the  two 

methods,  and  results  will  be  published  in  the  next  eight  months  (Ref.  1). 

Not«  Maniucript  submitUd  Dtccmbar  6,  1977. 
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As  a result  of  the  availability  of  computers  in  the  last  10  to  15  years  for 
developments  in  the  area  of  dynamic  response,  many  analytical  solutions  on 
the  dynamic  response  of  structures  have  evolved  based  on  assumed  mathematical 
models.  This  development  has  advanced  the  analysis  many  years  ahead  of  the 
knowledge  of  material  and  member  behavior  (which  serves  as  the  basis  for 
modeling),  thus  leaving  a broad  gap  between  the  two  areas.  Closure  of  this 
gap  is  essential  to  improve  and  to  check  finite  element  predictions. 

Eguipment  survival  is  of  paramount  importance  in  the  design  of 
U.S.  Navy  strategic  and  tactical  systems  against  weapon  effects,  and  fre- 
quently becomes  a driving  consideration  when  the  overall  cost  of  such  a 
system  is  considered.  As  is  to  be  expected,  every  effort  is  made  during  the 
conceptual  and  design  phases  to  minimize  system  equipment  costs  while  trying 
to  ensure  a high  probability  of  survival. 

The  expected  levels  of  shock  and  vibration  must  be  established 
before  deciding  how  to  protect  the  equipment.  Because  current  computer 
technology  limitations  cause  excessive  run  times  with  their  accompanying 
large  costs,  modeling  assumptions,  and  somewhat  awkward  structural 
and  b1 ast - load i ng  approximations,  reliable  predictions  of  the  shock  and 
vibration  environment  inside  protective  structures  are  limited.  The  equipment 
must  be  assessed  for  survivability  by  determining  failure,  malfunction,  and 
degradation  levels  (often  referred  to  as  "fragility")  in  terms  of  threat 
motions,  i.e.,  the  characteristic  acceleration  time  histories  expected,  rather 
than  simple  sine  waves,  saw-tooth  shock  pulses,  or  presently  used 
approximations  (Refs,  ?.  to  5)  • 

An  alternative  procedure  for  determining  the  internal  response  motions 
of  a protective  structure,  particularly  in  the  higher  frequencies,  is  to 
measure  the  impedances  of  the  structure  from  exterior  loading  surfaces  to 
internal  points  of  interest.  In  this  manner,  selected  points  on  the  struc- 
ture are  vibrated  with  a fixed-force  vibrator  as  the  frequency  of  vibration 
is  smoothly  changed  over  the  frequency  range  of  interest.  By  measuring  the 
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structural  response  at  selected  points  inside  and  on  the  structure,  a 
quantitative  measure  is  obtained  of  the  structure  response  to  the  vibratory 
forte.  The  complex  ratio  of  the  output  acceleration  to  the  input  force 
(magnitude  and  phase)  can  then  be  multiplied  by  an  arbitrary  input  force  to 
predict  sti'jctural  motions.  With  an  adequate  structural  description,  the 
motions  generated  at  equipment  mounting  by  an  attack  may  be  predicted,  and 
the  structural  path  can  be  analyzed  for  methods  mitigating  shock  transmission. 
When  the  appropriate  transmitted  shock  time  histories  at  the  mounting  points 
of  equipment  are  obtained,  the  equipment  can  be  hardened.  One  of  the  better 
and  more  practical  descriptions  of  the  structural  load  paths  is  impedance 
and  mobility.  The  severity  of  postulated  attacks  will  range  from  linear  to 
varying  degrees  of  nonlinear  response  of  a structure,  such  that  high-force- 
level  impedance  measurements  are  required.  The  high-force  levels  are  required 
for  impedance  measurements  in  a quasi-l inear  sense  or  for  the  formation  of 
nonlinear  functionals. 

Transient  shock  tests  on  equipment  and  systems  to  simulate  the 
motions  induced  by  a conventional  explosive  or  nuclear  attack  are  largely 
limited  to  single-axis  test  machines.  Further  limitations  exist  in  the  size 
and  weight  of  equipment  that  can  be  tested.  Simulating  multiaxis  loading  on 
large  equipment  with  many  degrees  of  freedom  represents  a difficult  problem, 
as  it  is  impractical  to  generate  continuously  varying  forces  of  sufficient 
magnitude.  This  problem  becomes  extremely  difficult  or  impossible  where 
in-place  or  field  tests  are  required.  On  the  other  hand,  short  duration 
forces  of  large  magnitudes  over  a wide  frequency  range  can  be  generated  by 
pulse  generators  (Refs.  6,  7).  Since  a discrete  number  of  pulses  superficially 
presents  an  appearance  quite  different  from  a continuous  input  excitation 
signal,  it  becomes  necessary  to  select  the  pulses  in  such  a way  that  the 
resulting  vibration  of  the  structure  matches  as  closely  as  possible  the 
response  (e.g.,  displacement,  velocity,  or  acceleration)  produced  by  the 
continuous  input  force,  as  determined  by  an  appropriate  error  criterion. 

Pulse  excitation  can  be  readily  demonstrated  analytically  to  induce  structural 
resf'onse  motions  closely  approximating  predicted  motions. 
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Numerous  physical  devices  can  be  applied  for  pulse  testing.  Their 
selection  and  use  depends  on  thrust  levels,  pulse  duration,  number  of  pulses, 
and  spatial  distribution  required.  The  class  of  pulse  generators  includes 
mechanical  (Ref.  8),  cold  gas,  gas/hydraulic,  hot  gas  (chemical),  and 
explosive  (point  charges).  Pulse  simulation  to  match  threat-criterion 
response  for  several  massive  weapons-systems  equipment  has  already  been 
accomplished  (Ref.  9).  Plans  for  pulse  testing  very  large  structural  systems 
to  damage/destruction  levels  are  now  being  made. 

The  combination  of  impedance  measurements,  response  predictions 
based  on  these  measurements,  and  pulse-simulation  testing  when  applied  to 
large  systems  offers  very  attractive  economic  alternatives  to  other  methods. 

For  example,  the  simulation  tests  of  Reference  9 provided  a satisfactory 
demonstration  of  adequacy  and  hardness  to  meet  mission  requirements  at  a 
cost  of  less  than  1%  of  what  was  expended  on  a different  weapon  system  for 
the  same  purpose.  In  this  other  weapon  system,  more  conventional  methods 
were  used,  and  the  results  were  less  than  conclusive.  Another  example  is 
the  present  construction  of  a $200,000,000,  two-axis,  50  ft  x 50  ft  shake- 
table  facility  by  the  Japanese  government  in  association  with  Japanese 
industry  to  test  nuclear-power  components  and  scale  models  to  earthquake 
levels  (Ref.  10).  Model  construction,  transport,  and  test  costs  would  be  an 
additional  expense. 

1.2  BACKGROUND 

1.2.1  UNITED  STATES  NAVY 

The  onset  of  the  development  of  nuclear-powered  submarines  generated 
a sustained  need  for  silencing.  One  of  the  principal  tools  for  the  analysis 
and  mitigation  of  structure-borne  sound  is  impedance.  Under  the  leadership 
of  the  Navy,  impedance  was  studied,  evaluated,  and  applied  by  naval  labora- 
tories, universities,  and  industry.  « 
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A.O.  Sykes  was  a major  contributor  at  the  Naval  Ships  Research 
and  Development  command,  as  was  Fred  Schloss  with  his  development  of  the 
impedance  head.  At  the  Naval  Research  Laboratory,  R.E.  Blake,  R.O.  Belsheim, 
G.J.  O'Hara,  I.  Vigness,  and  G.M.  Remmers  added  substantial  new  information, 
both  analytically  and  experimentally.  Through  the  Office  of  Naval  Research, 

N.  Perrone  provided  project  direction  on  impedance  in  both  industry  and  in 
the  academic  community. 

Analysis,  applications,  and  implementation  of  impedance  at  the 
Electric  Boat  Company  were  principally  carried  out  by  L.H.  Chen  and  I.P.  Vatz. 
R.  Plunkett  and  O.F.  Muster  conducted  similar  work  at  General  Electric.  The 
Pennsylvania  State  University  established  academic  leadership  in  impedance 
research,  particularly  in  its  associated  Ordnance  Research  Laboratory. 

Impressive  contributions  to  structural  dynamics  evolved  from  the 
above  activities  as  can  be  observed  from  the  following  structural  impedance 
I i St i ngs ; 


• Compiex  damping  and  loss  mechanisms 

• Th6venin  and  Norton  equivalents 

• Dynamic  absorbers 

• Isolation  with  finite  impedances 

• Component  connections 

• Lumped  parameter  systems 

• Distributed  systems 

• Beams  with  numerous  boundary  conditions 

• Built-up  structures 

• Plates  and  shel 1 s 

• Complex  and  statistical  force  fields 

• Systems  impedance  measurements  (ships  and  submarines) 
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The  listing  is  far  from  complete.  Published  work  originating 
from  this  list  is  extensive,  with  much  of  it  being  found  in  the  , Journal  of 
the  Aaoustiaal  Society  of  America  and  in  the  Shock  and  Vibration  Bulletin 
published  by  the  Shock  and  Vibration  Information  Center,  Naval  Research 
Laboratory . 


Fortunately,  many  of  the  fundamentals  of  U.S.  Navy  research  have 
been  preserved  and  consolidated  in  monographs  and  reference  books  such  as  the 
following:  Colloquium  on  Mechanical  Impedance  Methods,  R.  Plunkett,  ed . 

Amer.  Soc . of  Pech.  Engrs. , New  York,  Dec.  2,  1958;  J.C.  Snowden,  Vibration 
and  Shock  in  Damped  Mechanical  Systems,  John  Wiley  and  Sons,  New  York,  1988; 
E.  Skudrzyk,  Simple  and  Complex  Vibratory  Systems,  the  Pennsylvania  State 
University  Press,  University  Park,  PA,  1968;  V.H.  Neubert,  "Impedance  Related 
to  Modeling  and  Analysis  of  Structures,"  nnonograph.  Shock  and  Vibration 
Information  Center,  Naval  Research  Laboratory,  Washington,  D.C.  (to  be 
published  circa  1978). 

This  summary  of  the  impedance  activities  of  the  U.S.  Navy  omits 
many  contributors  only  because  of  the  need  for  brevity. 

1.2.2  EUROPEAN  WORK 

In  the  United  Kingdom,  impedance  work  on  structures  has  been  of  a 
long  and  consistent  duration,  commencing  after  its  formulation  by 
Professor  A.G.  Webster  of  Clark  University  in  191^.  Significant  work  was 
done  by  R.G.  White  and  F.  Kandianis  at  the  Institute  of  Sound  and  Vibration 
Research,  University  of  Southampton,  England,  where  time-series  digital  data 
analysis  was  extensively  applied  (Refs.  11  to  17).  The  work  at  Southampton 
was  sponsored  by  the  Royal  Navy.  Both  the  French  and  Swedish  Navies  are  also 
involved  in  the  application  of  impedance  techniques  for  submarines  (Ref.  18). 
Swedish  efforts  proceed  from  their  underwater  depth-charge  tests  on  a sub- 
marine section  (Ref.  19). 
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1.2.3  AEROSPACE,  MACHINE  TOOL,  AND  AUTOMOTIVE  INDUSTRIES 


The  introduction  of  the  Fast  Fourier  Transform  and  systems  identi- 
fication altered  the  direction  of  impedance  to  one  of  modal  analysis.  In 
this  approach,  data  measurements  are  processed  into  a form  compatible  with 
finite  element  modeling.  Numerous  d ig i ta 1 -process ing  methods  have  been 
developed  for  the  conversion  of  data  into  mode  shapes,  frequencies,  mass  and 
stiffness  matrices,  and  damping  representation  (Refs.  20  to  26).  Thus,  while 
an  impedance  measurement  is  a nonparamet r i c sample  of  a distributed  system, 
this  type  of  data  is  altered  to  the  form  of  uncoupled  normal  modes.  Tt v 
primary  objective  of  performing  modal  vibration  tests  is  to  develop,  verify, 
refine,  evaluate,  correlate  with  or  obtain  confidence  in,  the  analytical 
mathematical  model  and  the  dynamic  analyses. 

Applications  have  been  made  to  orbiting  vehicles,  missiles,  aircraft, 
ground  transportation,  nuclear  facilities,  automobiles,  and  power  plants. 

These  applications  have  been  facilitated  by  improvements  in  measurement  and 
testing  techniques,  sensor  developments,  and  on-line  digital  processing.  The 
latter  case  has  occurred  with  the  introduction  of  the  minicomputer.  These 
changes  have  yielded  improvement  in  measurement  accuracy,  and  a very  substan- 
tial reduction  in  costs. 

1.2. A DEFENSE  NUCLEAR  AGENCY  AND  U.S.  ARMY  CORPS  OF  ENGINEERS 

The  need  for  accurate  environmental/structural  motions  over  a wide- 
frequency  band,  caused  by  nuclear  weapon  effects,  motivated  C.C.  Huang,  Corps 
of  Engineers,  Huntsville,  and  J.G.  Lewis  and  E.P.  Sevin  of  the  Defense  Nuclear 
Agency  to  direct  projects  where  structural  responses  could  be  computed  from 
impedance  measurements.  The  Corps  of  Engineers'  projects  involved  the  motion 
prediction  of  massive  equipment  (200,000  lb)  from  input-mounting  locations 
to  cr i t ica I - i n terna 1 locations  of  the  equipment  as  determined  by  transfer 
functions.  The  Defense  Nuclear  Agency  projects  involved  motion  predictions 
of  aboveground  and  flush-buried  protective  structures  as  determined  by 
impedance  functions.  These  projects  are  summarized  in  Section  2. 
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The  Corps  of  Engineers,  Huntsville,  also  required  in-place  testing 
of  large  massive  equipment  where  the  test  motions  would  reasonably  match 
the  acceleration  time  histories  predicted  from  transfer  functions.  To  meet 
this  need,  a mechan ica 1 -phase  generator  was  developed  and  successfully 
applied.  Details  of  this  pulse  system  were  briefly  discussed  in  Section  1.1, 
and  are  elaborated  in  Section  2. 
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SECTION  2 


RESPONSE  PREDICTIONS  AND  INDUCED  MOTION  SIMULATION 

FROM 

MEASURED  AND  SCALED  SYSTEM  FUNCTIONS' 

2.1  SYSTEM  FUNCTIONS 

2.1.1  DEFINITIONS 

Dynamic  properties  of  passive  systems  may  be  represented  by  a 
complex  ratio  of  the  input  loading  to  the  output  response  of  a system.  This 
complex  ratio  may  be  presented  in  the  form  of  magnitude  and  phase,  real  and 
imaginary,  and  real  vs.  imaginary,  all  as  a function  of  frequency.  Inverse 
transformation  provides  yet  another  form  of  this  ratio  in  the  time  domain, 
the  impulse  function.  The  compiex  ratios  obtained  may  be  viewed  as  samples 
from  a continuum  both  in  space  and  bandwidth,  and  are  nonparametric.  These 
functions  can  be  determined  analytically  from  a mathematical  model  or  can  be 
physically  measured.  Figure  1 illustrates  the  concept: 


I nput 
O— 


Structure 


Output 


• Impedance  as  function  of  frequency 

7(  ) m Input  Force 

Output  Ve  loci  ty 

• Transfer  function  as  function  of  frequency 

) a Output  Motion  ^ Output  Force 
' Input  Motion  Input  Force 


FIGURE  I.  IMPEDANCE  AND  TRANSFER  FUNCTIONS  CONCEPTUALIZED 
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Mechanical  impedance  is  an  indicator  of  how  a structure  responds 
to  a vibratory  force.  The  response  of  a vibrated  point  on  a structure  is 
inversely  proportional  to  the  impedance,  which  will  be  defined  herein  for  a 
given  frequency  as  the  driving  forre  on  the  structure  divided  by  the  peak 
velocity  at  a point  on  the  structure.  The  motion  of  the  structure  may  be 
recorded  as  displacement,  velocity,  or  ac'e I era t i on , and,  when  ratioed  with 
input  force,  the  various  system  functions  are  defined  by  the  following  terms: 


Force 
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I nertance 


Another  important  function  is  the  ratio  of  input  motion  to  output 
motion  or  of  input  force  to  output  force.  This  is  defined  as  a transfer 
function.  Thus,  if  the  base  of  a structure  is  put  into  motion  by  a blast  or 
earthquake,  multiplying  this  base  motion  by  the  transfer  function  produces 
the  motion  at  another  part  of  the  structure.  The  terms  "impedance"  and 
"transfer  function"  have  been  used  to  simplify  the  presentation  of  this 
report;  and  where  specific  functions  are  displayed,  as  in  data  plots,  these 
terms  are  defined. 


2.1.2  EFFECTS  OF  NOISE  AND  SIGNAL  PROCESSING 

The  total  output  of  the  data  measurement  device  is  the  signal. 

This  output  signal  passes  into  a recording  system  that  produces  i signature 
on  magnetic  tape.  Subsequently,  the  analog  signature  is  converted  to  a 
digital  signature  also  on  magnetic  tape.  The  digitt-i  data  are  then  in  a 
form  suitable  for  analysis  by  digital  processors.  In  this  chain  of  events,  • 

however,  there  are  a number  of  opportunities  for  electronic  and  mechanical 
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errors  to  occur.  These  errors  are  commonly  defined  as  noise  or  drift  and 
offset.  Data  measurements  to  obtain  transfer  and  impedance  functions  are 
very  susceptible  to  noise,  in  both  acquisition  and  processing  operations. 

These  spurious  effects  are  often  of  a magnitude  sufficient  to  highly  distort 
the  measured  signal.  A display  of  signal  and  noise  (Ref.  27)  is  shown  in 
Figure  2 for  a sine  wave  input  signal  sweeping  from  1*0  to  80  Hz. 

Complex  ratios  of  signals  containing  moderate  amounts  of  noise 
exhibit  a pattern  similar  to  Figure  3 for  the  inertance  magnitude  function 
lX/F(ju))].  Noise  both  enhances  and  reduces  resonant  peaks  with  no  clear 
pattern.  A more  sensitive  indication  of  noise  can  be.  found  in  phase  plots, 
i(i(juj),  of  Figure  k when  the  phase  is  presented  for  0 to  -nn  deg  rather  than 
in  the  conventional  form  of  tl80  deg.  For  noise-free,  mu  1 1 i resonant  linear 
systems,  the  phase  plot  should  be  essentially  a monotonic  line  of  phase  as  a 
function  of  frequency.  The  slope  of  this  curve  provides  the  arrival  time  of 
the  signal,  which  is  given  as: 

T - i 

iii 

whe  re 

1 “ Arrival  time  (sec) 

<p  “ Phase  change  (radians) 
u)  = 27if  frequency  (radians) 

Examination  of  Figure  ^a  shows  an  erratic  pattern  of  phase,  and 
this  phase  pattern  can  be  observed  in  the  plot  of  the  impulse  function  of 
Figure  5a,  The  phase  plot  of  Figure  4b  tends  to  be  a better  signal,  yet  the 
impulse  function  of  Figure  5b  exhibits  a high  degree  of  noise  and  some 
effects  of  symmetry  at  the  end  of  the  trace.  Obviously,  if  these  data  in 
Figures  A and  5 were  to  be  used  in  their  present  forms,  very  poor  computations 
• would  result. 


(Text  continued  on  p.  16.) 
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FIGURE  2.  DISPLAYS  OF  SIGNAL  AND  NOISE  FOR  DATA  ANALYSIS 
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FIGURE  3.  example  OF  INERTANCE  MAGNITUDE  AFFECTED  BY  NOISE 
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FIGURE 


(a)  Transfer  function  phase  of  equipment 


(b)  Inertance  function  phase  of  a protective  structure 
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Another  method  for  the  detection  of  noise  is  the  coherence  function 
(Ref.  28),  which  is  expressed  by  the  following: 


|G  (f)l' 

nffxTTT  : ' 

y 


where 


G (f) 

X 


G (f) 
xy 


G^(f) 


Power  spectral  density  function,  Fourier  transform  of 
autocorreiation  function  of  x(t) 


Cross-spectral  density  function,  Fourier  transform  of 
cross-correlation  function  of  x(t)  and  y(t) 

Power  spectral  density  function,  Fourier  transform  of 
autocorrelation  function  of  y(t) 

Coherence  function 


2 

When  > (f) 

xy 

said  to  be  incoherent 
then  x(t)  and  y(t) 


“0  at  a part 
(uncorrelated) . 
are  said  to  be 


icular  frequency. 

When  y^  (f)  = 
'xy 

fully  coherent . 


x{t)  and  y(t)  are 
1 for  all  frequencies. 

If  the  coherence  function 


is  greater  than  zero  but  less  than  unity,  one  of  three  possible  conditions 
exist : 


• Extraneous  noise  is  present. 

• The  system  relating  x(t)  and  y(t)  is  not  linear. 

• y(t)  is  an  output  due  to  an  input  x(t)  and  other 
inputs. 

Extraneous  noise  suppresses  the  coherence  fucntion,  and  the  display 
of  Figure  6 serves  as  a detector  for  noise. 
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2.1.3  REDUCTION  OF  NOISE  AND  ERROR  SIGNALS 


2. 1.3-1  I nstrumentat ion 

It  goes  without  saying  that  adequate  and  proper  instrumentation  is 
essential  to  measure  and  record  data  with  minimal  noise  and  error  signals. 
This  involves  systematic  methods  of  grounding,  shielding,  and  bonding  for 
electromagnetic  compatibility  and  interference.  Calibration  of  the  system 
in  bandwidth  and  amplitude  is  required  prior  to  each  test  run  and  on-line 
monitoring  of  signals  is  required  during  tests.  All  transducers  must  be 
checked  for  cross-coupling  and  acoustic  pickup.  Experienced  professionals 
with  a proven  performance  record  in  instrumentation  must  be  selected. 

2. 1.3.2  Signal  to  Noise 

Input-force  signal  levels  should  be  set  at  the  highest  level  com- 
patible with  test  objectives  and  the  linearity  of  the  test  article  to  obtain 
more  favorable  signal-to-noi se  ratios.  This  objective  is  more  often  limited 
to  the  force-generating  equipment  available.  Impedance  tests  performed  on  a 
1/12-scale  model  of  a protective  structure  at  the  Defence  Research  Establish- 
ment, Suffield,  Canada,  ranged  from  15  to  30  lb,  while  identical  impedance 
measurements  on  the  prototype  structure  in  North  Dakota  ranged  from  1000  to 
5000  1b.  Nominally,  this  change  in  force  level  between  the  two  structures 
corresponded  to  the  scaled  impedance  difference  of  lAA. 

Checks  can  be  made  to  ensure  that  the  system  ambient  electrical 
and  mechanical  noise  levels  do  not  interfere  with  the  impedance  data.  This 
is  accomplished  by  making  plots  of  impedance  vs.  frequency,  using  as  data  the 
ambient  signal  level  from  an  accelerometer  mounted  on  the  structure  and  a 
force  signal  level  from  the  operating  but  mechanically  isolated  vibrator. 

The  vibrator  can  be  isolated  by  suspending  it  on  cables  and  attaching  a 
reaction  mass;  or  if  the  size  of  the  vibrators  used  preclude  such  direct 
isolation,  than  the  force  signal  can  be  electrically  simulated.  Figures  7 
and  8 for  a prototype  protective  structure  and  its  1/12-sca1e  model  show 
noise  overlaid  on  impedance  data. 
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FIGURE  7.  TYPICAL  PROTOTYPE  PROTECTIVE  STRUCTURE  TRANSFER 
IMPEDANCE/NOISE  COMPARISON 


19 


MECHANICAL  IMPEDANCE,  LB-SEC/IM. 


FREQUENCY.  HZ 

FIGURE  8.  1/12-SCALE  MODEL  PROTECTIVE  STRUCTURE  I MPEOANCE/NO I SE  COMPARISON 
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2. 1.3-3  Phase-Shifting  Errors 

With  magnetic  tape  recording  machines,  phase  shifts  occur  even  though 
the  heads  are  within  manufacturer's  tolerance.  Recording  data  on  one  machine 
and  playing  back  on  another  compounds  the  problem.  Phase  shifts  vary  directly 
with  frequency  and  inversely  with  tape  speed  as  per  the  following  expression: 

Head 

Phase  shift  error  = — * ' gnmen  t , i n ._  ^ frequency  x 360  deg 

Tape  Speed,  i ps 

Tape  machines  should  be  calibrated  for  each  channel  with  an 
oscillator  over  the  frequency  band  of  interest.  Phase  corrections  may  be 
used  subsequently  when  the  data  have  been  digitized. 

2.1.3.^  Drift  and  Offset  ^ 

Offsets  and  drifts  occur  in  many  cases  and  are  independent  of  the 
measurement  instrument  response  to  externally  applied  loads.  In  such  cases, 
record  offsets  or  drifts  are  stationary,  i.e.,  the  same  error  is  present 
before,  during,  and  after  the  transient  era.  A simple  means  of  correcting 
these  errors  is  to  determine  the  offset  and/or  drift  on  a time  segment  before 
the  shock  arrival  and  to  extrapolate  the  appropriate  correction  to  the 
transient  segment. 

Experience  has  shown  that  late-time  detrending  is  the  most  reliable 
way  of  correcting  most  baseline  errors.  The  technique  (Ref.  29)  is  illustrated 
in  Figure  9-  In  many  cases,  a linear  least-squares  fit  to  the  late-time  data 
is  adequate  to  achieving  the  desired  correction.  When  the  record  length  is 
toe  short  to  use  the  preceding  method,  the  baseline  error  is  corrected  by 
the  early-time  procedure  illustrated  in  Figure  10.  Again,  linear  detrending 
is  adequate  to  removing  the  baseline's  errors.  The  line  fitted  to  the 
early- time  segment  is  defined  as 

q ( t)  = a + bt 


Occasionally,  the  data  justify  a higher  order  correction,  such  as  a parabolic 
fit  to  the  data . 
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FIGURE  10.  EARLY  TIME  DETRENDING  OF  DATA 
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where  a is  the  intercept  and  b is  the  slope.  The  constants  a and  b 
are  evaluated  by  mi  n i m i ^ i ncj  the  function 

Q • jf(l)-9(t)|^ 

such  that 


where  f(t)  is  the  data  to  be  fitted.  It  Is  to  be  noted  that  the  first 
equation  considers  both  data  offset  and  drift.  If  only  offset  is  to  be 
removed,  the  constant  b vanishes.  Tests  for  effectiveness  of  the  offset 
and  detrending  correction  is  by  integration  to  determine  residual  velocity 
and  displacement  and  compare  with  observed  or  known  results. 

2 . 1 . 3 • 5 Discrete  Noise  Rejection 

The  presence  of  discrete  noise  is  identified  from  the  Fourier 
transform  of  the  record  where  no  signal  Is  present.  Figure  11  shows  an 
example  of  such  a transform.  The  fundamental  60  Hz  power-supply  noise  and 
the  even  harmonics  out  to  the  8th  harmonic  are  clearly  evident.  Also  noted 
in  Figure  II  is  the  random  component  of  noise  exemplified  by  the  nearly 
constant  gain  across  the  entire  frequency  band. 

The  removal  of  discrete  noise  can  be  accomplished  in  various  ways, 
notably  by  band-reject  filtering  and  time-domain  subtraction.  The  former 
technique  is  tedious  and  unreliable  because  band-reject  filtering  can  remove 
important  signal  components.  The  latter  technique  is  potentially  more 
accurate  as  described  below. 
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FIGURE  II.  ILLUSTRATION  OF  DISCRETE  AND  RANDOM  NOISE 


The  time-domain  subtraction  technique  (Ref.  29)  consists  of  sub- 
tracting the  discrete  noise  segment  of  a record  from  the  s i gna I -p 1 us-noi se 
segment.  Two  major  changes  will  result  from  this  process.  First,  if  the 
noise  is  properly  phased  with  the  signal  plus  noise,  the  discrete  noise 
components  will  cancel.  Implicit  in  this  result  is  the  assumption  and  fact 
that  the  entire  measurement  and  playback  systems  are  linear  for  the  measure- 
ment voltages.  Second,  the  random  component  of  noise  will  be  amplified 
because  the  subtraction  of  two  random  segments  will  result  in  a record  whose 

variance  is  larger  than  for  either  of  the  segments.  The  penalty  of  the 

second  result  is  accepted  in  the  interest  of  achieving  the  benefits  of  the 
first  resu I t . 

To  establish  the  proper  phase  times  between  the  noise  and  signal- 
plus-noise  segments  of  a record,  a procedure  of  band-pass  filtering  and 
cross  correlation  is  implemented.  Band-pass  filtering  at  the  fundamental 
discrete  noise  component  (60  Hz)  and  at  the  highest  discernible  harmonic 
thereof  is  performed  on  an  entire  record.  This  procedure  eliminates  all  fre- 
quencies other  than  those  in  the  immediate  neighborhood  of  the  fundamental 
and  the  harmonic.  A noise  portion  of  the  record  is  subsequently  cross 

correlated  with  the  signal  portion  of  the  record  via  the  digital  Implementa- 

tion of  the  equation  below 


R (t) 

XX 

n s+n 


I im  _1_ 
T~-  T 


f X (t)  X (t  + t)  dt 
Jo  n s+n 


where  x^  and  denote  the  band-passed  noise  and  signal  plus  noise, 

respectively.  The  cross-correlation  function  derived  from  evaluating  the 

equation  for  successive  lag  times  x ultimately  establishes  the  particular 

delay  when  both  the  fundamental  and  harmonic  are  in-phase  between  x and 

n 

Taking  into  account  the  particular  i thus  derived,  the  unfiltered 
noise  segment  of  the  record  can  be  subtracted  from  the  unfiltered  signal-plus- 
noise  segment  to  achieve  the  desired  result,  i.e.,  a record  purged  of  the 
discrete  noise. 
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2. 1.3.6  Kandom  Noise 


Random  noise  of  the  same  bandwidth  as  the  signal  presents  special 
problems,  particularly  for  single  tests.  When  repetitive  tests  can  be  made, 
averaging  will  substantially  reduce  the  random  noise  content. 

Repetitive  tests  of  the  same  measurement  can  be  made,  and  transfer 
or  impedance  functions  are  formed  for  each  test.  These  functions  are  treated 
as  repetitive  waveforms  and  averaged.  The  variance  of  the  noise  present  in  the 
function  is  reduced  by  the  number  of  samples  and  consequently  the  rms  value 
of  the  noise  is  reduced  by  the  square  root  of  the  number  of  samples  averaged. 
This  method  assumes  that  the  mean  value  of  the  noise  is  zero,  that  the  noise 
is  a random  variable,  and  that  each  sample  (noise)  is  statistically  independent 
An  example  of  this  procedure  is  shown  in  Figure  12  for  20  averages.  This 
figure  should  be  compared  to  Figure  3.  which  represents  a single  measurement 
of  the  same  test  article. 

Repetitive  system  function  averaging  as  well  as  single  measurement 
tests  can  be  supplemented  by  curve  fitting  to  the  data.  This  curve-fitting 
method  is  a very  useful  procedure  and  produces  response  predictions  that  are 
often  reasonable  and  not  excessively  in  error  even  when  used  with  system 
functions  containing  noise.  In  this  latter  case,  technical  judgments  and 
end  objectives  of  the  project  must  determine  whether  to  use  curve  fitting. 

Other  noise  reduction  and  evaluation  methods  have  been  developed 
and  are  useful  in  varying  degrees  for  single  test  measurements.  These  include 
PSD  (Ref.  29),  Fourier  transform  (Ref.  29),  phase  optimization  methods 
(Ref.  30),  and  coherence  (Ref.  31).  At  this  time,  these  methods  are 
approximate  and  costly  to  implement. 

2.1.1*  excitation  OF  STRUCTURES 

Selection  of  driving  functions  are  constrained  by  the  physical  hardwa 
available  to  generate  the  required  force  levels,  shaping  characteristics,  and 
bandwidth.  The  most  commonly  used  force-function  generators  are  electrodynamic 


and  electrohydraul ic  shakers  that  produce  continuous  periodic  (usually 
sinusoidal)  functions  that  are  controllable  in  amplitude,  time  duration, 
and  frequency.  With  sophisticated  control  systems,  these  devices  are  also 
capable  of  generating  random  and  complex  force  functions  as  well  as  simple 
pulses.  Rotating  eccentric  mass  shakers  are  also  used,  primarily  for  civil 
structures. 

Shock  machines,  impact  hammers,  pulse  generators, 
produce  a second  family  of  driving  functions,  the  transient 
types  of  functions  may  include  many  pulses  of  varying  forms 
triangle,  square),  depending  on  the  required  input. 

Methods  of  excitation  that  have  come  into  common  use  and  are 
discussed  in  this  section  are 

• Slow  sine  sweep 

• Random 

• Rapid  sine  sweep  (chirp) 

• Impulse 

• Pulse  train 

2.I.I4.I  Slow  Sine  Sweep 

In  the  past,  the  frequency  response  characteristics  of  structures 
have  been  measured  by  quasi -steady-state  methods.  One  method,  referred  to 
as  the  slow  sine-sweep  method,  involves  varying  an  input  sinusoidal  forcing 
function  over  the  frequency  region  of  interest  and  measuring  the  output 
responses.  It  has  been  used  extensively  in  impedance  and  transfer  function 
testing  (Refs.  30  and  32).  Considerable  error  may  occur  if  steady-state 
response  is  arbitrarily  assumed.  To  minimize  such  error,  convergence  tech- 
niques are  used  in  the  testing  procedures.  Convergence  techniques  involve 
progressively  slower  sweep  rates  (or  t ime-of-f requency  variation)  until 
essentially  asymptotic  response  amplitudes  are  obtained. 


and  explosives 
functions.  These 
(e . g . , half-sine. 
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A/D  conversion  and  the  subsequent  computer  processing  of  slow  sine- 
sweep  data  generate  cost  and  technical  problems  owing  to  the  large  amount  of 
data  required.  To  circumvent  this  problem,  sophisticated  A/D  conversion 
techniques  have  been  developed  that  allow  for  selective  windowing  of  the  real- 
time response  data  (Refs.  30  and  32).  The  window-sampled  data  are  processed 
for  frequency,  magnitude,  and  phase,  and  these  data  are  stored  as  sets  of 
digital  channel  elements.  For  many  large  data  acquisition  projects,  however, 
this  processing  method  is  impractical  due  to  the  large  time  and  cost  require- 
ments. In  addition,  averaging  the  system  functions  to  reduce  noise  is  not 
practical . 

The  slow  sine  excitation  is  useful  for  quick-look  field  data 
processing  with  analog  impedance  equipment.  On  large  structures,  this 
method  permits  resonance  buildup  of  the  structure,  thereby  generating  more 
favorable  s i gna 1 - to-noi se  ratios.  Vibration  machines  ranging  from  a few 
pounds  to  50,000  lb  of  output  force  are  available.  Both  low  and  high  fre- 
quency rolloff  of  force  occurs  and  is  a function  of  the  design  characteristics 
of  each  vibrator. 

2.1.4. 2 Random  Function 

The  random  driving  function  can  include  all  frequencies  of  interest 
over  a selected  range.  It  is  defined  as  an  uncorrelated  function  for  which 
the  mean  is  zero  and  the  power  spectral  density  (PSD)  is  constant  or  shaped 
over  the  frequency  region  of  interest. 

'*5°  ■ ^ ■ b'o  T " f 

in  which 

B “ Frequency  bandwidth 

“ Power  spectral  density 

a “ rms  of  random-force  input 
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System  functions  (Ref.  28)  may  be  measured  from  input-output  records 
of  random  excitation  by  the  following: 


-0  (f)  ■.lf\ 

|G.  {f)|e  - |H(f)|e'-'^^^^  G (f) 


xy 


lH(f) 


|G  (f) 
xy 

G^(f) 


' <'(f) 


where 


Gx(f) 

'^xy<^) 

0,  (f) 

xy 

♦ (f) 
lH(f)| 


Input  spectral  density 
Cross-spectral  density  (input/output) 
Cross-spectral  density  phase  factor 
Phase  system  function 
System  function  magnitude 


This  type  of  excitation  permits  sampling  of  the  data  records  to  form 
a succession  of  system  functions,  which  when  averaged  substantially  reduce 
noise.  Figure  12  is  an  averaged  system  function  from  20  trials  using  random 
excitation.  Input  power  per  frequency  from  the  vibrator  is  reduced  for 
broad  band  random  compared  to  the  slow  sine  excitation  method.  Quick- 
look  at  the  system  function  must  be  delayed  until  the  records  have  been 
processed,  unless  a minicomputer  processor  is  available  (Time-Data,  HP, 
Spectral  Dynamics).  Force  input  rolloffs  at  both  low  and  high  frequencies 
are  functions  of  the  particular  vibrator.  Figure  13a  is  an  example  of  a 
random  force-time  history. 
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2. 1.^.3  Rapid  Sine  Sweep  (CHIRP) 


Generally  similar  to  the  slow  sine-sweep  function  is  the  rapid  sine- 
sweep  (chirp)  function,  which  differs  mainly  in  the  time  duration  of  frequency 


sweeping.  Figure 

13b  is  a 

force-time  history  upsweep  chirp.  Two  forms  of 

this  f unc  t i on 

are 

common  1 y 

used  in  testing,  one  from  Reference  33. 

F(t) 

s 

2 

Fo  sin  (at  + bt) , 0 < t < T,  increasing  frequency  sweep 

where 

Fo 

s 

Maxi  mum 

force  amplitude 

a 

= 

■ '2  > 

T 

Total  record  length,  sec 

^ = 

Initial  frequency,  Hz 

^ * 

Final  frequency,  Hz 

b 

2n  f. 

and  the  second  from  Reference  3^: 


F(t) 


Fo  sin 


-2"  N-  log  (l  - 


frequency  sweep 


and 


, for  increasing 


F(t) 


Fo  s i n 


2n  N'  log^  (l  + M) 


frequency  sweep, 


, for  decrees i 


ng 


where 

Fo  “ Maximum  force  amplitude 

N'  ■ Effective  number  of  cycles  at  each  frequency 
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fo  • Initial  frequency,  Hz 
f ” Final  frequency,  Hz 
t “ Sweep  time  from  fo  to  f 

» N'  (-^  - y)  , increasing  frequency 
“ "N'  - j)  , decreasing  frequency 

Chirp  testing  may  be  performed  in  either  direction,  from  low  to  high  frequency 
or  the  reverse.  As  in  random  testing,  short-record  lengths  may  be  obtained 
for  the  frequency  bandwidth  of  interest  with  consequent  ease  in  digitizing 
the  record.  Function  generators  controlling  the  vibrators  may  be  programmed 
for  numerous  chirp  sweeps  to  permit  averaging  of  the  system  functions  for 
noise  reduction.  For  chirp  testing,  additional  control  systems  may  be  required 
for  the  vibrator  to  maintain  programmed  input.  An  example  of  a single  chirp 
sweep  is  the  system  function  shown  in  Figure  3-  Limitations  of  quick-look 
data  and  frequency  rolloff  of  the  vibrator  are  the  same  as  for  random 
excitation.  Major  advantage  of  the  chirp  test  is  the  higher  force  input 
attainable  per  frequency. 

2 

Chirp  testing  in  the  form  of  F^  sin  (at  +bt)  provides  a reasonably 
constant  force  input  spectrum  as  may  be  observed  in  Figure  Ilf.  For  general 
testing  this  is  the  desired  form.  The  log  sweep  form  provides  the  frequency 
spectrum  given  in  Figure  15.  In  this  log  sweep  form,  the  low  spectrum 
amplitudes  at  the  high  frequency  yield  undesirable  s i gna I -to-noi se  ratios. 

The  log  sweep  chirp  is  desired  by  some  investigators  for  the  form  or  shape 
of  the  shock  sp'^ctrum,  which  may  be  generated. 

The  shock  .pe  rrum  is  an  envelope  of  the  response  amplitudes  of  n 
single-degree-of-freedum  (SDOF)  oscillators  as  n * The  amplitudes  are 
dependent  on  the  nature  of  the  input  and  the  damping  assumed  for  the  oscilla- 
tors. The  shock  spectrum  is  used  to  determine  how  the  oscillators  (and  the 
modes  for  real  structures)  are  excited  towards  their  maximum  amplitudes.  It 
is  not  enough  just  to  have  a frequency  present  in  the  input.  It  must  be 


TUOE,  IN  THE  FORM  OF 


o 

o 
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FIGURE  15.  CHIRP  TEST  FOURIER  MAGNITUDE  IN  THE  FORM  OF  A LOG  SWEEP 


present  with  adequate  osciiiations  to  fully  excite  the  mode.  The  evaluation 
of  whether  or  not  an  input  will  produce  this  desired  effect  is  accomplished 
by  comparing  the  shock  spectrum  for  several  values  of  damping.  As  Figure  16 
shows,  the  number  of  cycles  required  is  heavily  dependent  on  the  value  of 
Q,  (C/Cc  “ t “ 1/2Q).  For  Q « 5(c  “ 0.10),  the  response  amplitude  has 
essentially  reached  a maximum  of  5 cycles,  but  for  Q • 25  (r,  » 0.02),  the 
maximum  has  not  been  reached  at  20  cycles.  Q ■ 5 is  a reasonable  value  to 
assume  for  structural  modes  and,  therefore,  from  Figures  16  and  17  it  can  be 
concluded  that  whenever  > 2,  the  input  is  adequate  to  fully  excite 

the  structural  modes. 

Shock  spectra  calculated  at  values  of  Q,  ” 5 end  Q “ 25  for  the 
chirp  (log  sweep)  driving  function  is  shown  in  Figure  18  from  which  the 
ratio  Q2^/Qg  will  give  an  effective  number  of  osciiiations  of  15-  This 
shock  spectrum  technique  is  a useful  tool  in  selecting  sweep  times  for 
chirp  testing. 


2.1. A. A Impulse  Functions 

The  impulse-function  family  comprises  a number  of  types.  Two  of 
the  types  that  are  mechanically  producible  are  the  short  duration,  single- 
impulse functions  from  shock  machines,  and  periodic  impulses  from  impact 
hammers.  The  simplicity  and  portability  of  some  of  the  devices  used  to 
generate  impulses  make  them  advantageous  for  many  engineering  applications. 
A typical  impulse  function  and  associated  Fourier  transform  magnitude  are 
presented  in  Figures  13c  and  19. 


A single  Impulse  function  that  can  populate  a high-frequency  region 

must  be  of  short  duration  as  seen  in  Figure  19  (e.g.,  let  F ■ 1/2  — » 

’ max  a 

500  Hz;  then  a ^ 5 ms).  The  short  duration  imposes  practical  restrictions 
on  shock  devices  in  formulating  acceptable  pulses.  For  lightweight  test 
articles,  a shock  hammer  was  developed  at  the  Naval  Research  Laboratories 
(Ref.  35)  circa  1958-  Practical  use  of  this  device  did  not  come  about  until 
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(Ref.  3'*) 


FIGURE  18.  SHOCK  SPECTRA  FOR  LOG  SWEEP  CHIRP  AT  Q * 25  AND 


FIGURE  19 
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ifLLiill/  vvilh  ttii.-  .idvi'nt  of  digital  spectrum  analyzers.  Today,  one  can  flail 
ui.us  (i'  iiu,  tiuiiKMOus  repeal  tests  to  average  down  noise  inherent  in  this 
e . i til  Ion  thod.  Given  sufficient  weigiit  or  force,  impulse  excitation  can 
he  .1  ..  .1  lo  , KL  i te  and  measure  very  low  frequencies. 

2 . 1 . ft . b fu  I se  Traots 

Specially  designed  mechanical  force  pulse  generators  have  been 
used  quite  Micces s f u 1 I y to  simulate  known  acceleration  excitations  for  vibra- 
tion testing  (Kef.  6).  These  generators  produce  pulse  trains  with  specified 
amplitudes  and  time  durations.  The  pulse  trains  can  be  optimized  via  computer 
stWtwaie  to  simulate  a specific  output.  The  pulse  trains  are  designed  to 
produce  a predetermined  Fourier  spectrum  shape  throughout  the  frequency  region 
of  interest  (Ref.  36).  The  pulse-train  method  is  particularly  advantageous 
because  of  the  portability  of  test  equipment,  the  short  duration  of  the 
inputs,  and  the  large  force  input  capabilities.  A typical  pulse  train  is 
presented  in  Figure  13d,  and  its  Fourier  magnitude  is  presented  in  Figure  20. 

Tfie  impulse  method  is  extended  by  the  pulse  train  to  achieve  greater 
nxit.  i I i zdt  ion  ot  the  test  article  by  pulse  sequencing  and  by  spectrum  shaping. 
Addi t ional I y,  a reasonable  degree  of  input  force  control  is  obtainable.  Force 
Ic'fls  of  20,000  lb  have  been  obtained,  and  100,000-lb  inputs  appear  achiev- 
able with  pi esent  devices.  The  pulse  train  method  of  excitation  provides  a 
high  force,  low  frequency  excitation,  a moderate  bandwidth,  and  a reasonable 
numbei  of  oscillations  per  frequency.  Quick- look  evaluation  is  restricted 
to  oscillograph  time  histories.  System  functions  must  be  processed  digitally. 

2.1.5  SYSTtM  SURVEY 

At  the  onset  of  a test,  a survey  is  recommended  to  ascertain  the 
linearity  of  the  system  by  multilevel  tests  and  reciprocity,  system  symmetry, 
and  cross-axes  coupling.  Quantitative  information  gained  at  an  early  stage 


FIGURE  20.  FREQUENCY  TRANSFORM  OF  FORCE-PULSF  INPUT 


can  permit  reductions  in  test  efforts  or  provide  for  alternative  procedures. 
Quick- look  on-line  plots  of  system  functions  are  a necessity. 


2. 1.5. 1 Multi  level  Tests 

Linear  systems  will  provide  system  functions  that  are  virtuai  over- 
lays of  each  other,  although  made  at  several  different  input  force  levels.  For 
the  Perimeter  Acquisition  Radar  Building  (PARB)  of  the  SAFEGUARD  system,  force 
input  of  1000  Ib  to  7000  Ib  demonstrated  for  this  range  that  essential 
linearity  of  the  building  existed.  Transfer  function  tests  across  isolators 
supporting  a 200,000-lb  control  room  platform  showed  nonlinearity  existing, 
as  may  be  observed  in  Figure  21.  Another  form  of  nonlinearity  (Fig.  22)  was 
uncovered  for  shock- i sol ated  platform  from  static  measurements  and  measure- 
ments made  while  the  platform  was  oscillated  at  1 Hz.  Subsequently,  response 
predictions  using  the  higher  transmission  values  indicated  that  the  local 
environment  was  well  below  the  hardness  value  of  the  nearby  mounted  equipment. 

Nonlinearity  measurements  of  transfer  inertances  along  a missile- 
silo  wall  sited  in  backfilled  soil  are  shown  in  Figure  23.  For  the  most  part, 
they  show  an  increase  in  compliance  with  increased  test  load. 

2. 1.5. 2 Rec iproc i ty 

Reciprocity  measurements  were  made  on  the  PARB  model  in  Canada  and 
the  prototype  PARB  in  North  Oal.ota  in  accordance  with  the  diagrams  in  Figure  2k. 
Reciprocal  impedance  measurements  for  both  buildings  are  given  in  Figures  25 
and  26.  Extensive  use  of  reciprocity  was  made  in  the  prototype  building,  as 
It  was  much  easier  to  move  accelerometers  around  the  exterior  surfaces  than 
the  large  and  heavy  vibrator. 

Measurements  on  a moderate  size  rei nforced-concrete  protective 
structure  buried  flush  to  ground  surface  disclosed  significant  noni i near i t i es 
(Fig.  27). 


(Text  continued  on  p.  51 •) 
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FIGURE  21.  PLATFORM  PARPP-CR:  TRANSFER  FUNCTION  MODULUS  (ANALOG)  FORCE  RATIO  ACROSS 
SPRING  ISOLATOR  AT  INPUT  FORCE  LEVELS  OF  200  LB  AND  600  LB 


FIGURE  22.  PLATFORM  MSRB-I25B:  TRANSFER  FUNCTION  MODULUS  ACROSS  PNEUMATIC  ISOLATOR, 
COMPARISON  OF  MEASUREMENTS  WITH  AND  WITHOUT  1-HZ  CONSTANT  OSCILLATION 


FIGURE  23.  TRANSFER  INERTANCE,  X/F(jw) , MAGNITUDE  FUNCTION  ALONG  A Ml  SSI LE 
SILO  WALL  SITED  IN  BACKFILLED  AND  COMPACTED  SOIL 


(a)  Motfel  ?AS3:  OIrecc  Tieasu remen :s  (frcm  surface  to  ■'espcrse  oc 


FIGURE  21*.  MEASUREMENT  METHODS  USED  ON  PARS  MODEL  AND  PROTOTYPE 


FIGURE  25.  MODEL  PARB  DATA  FOR  DETERMINING  FIGURE  26.  PROTOTYPE  PARS  DATA  FOR  DETERMINING 

RECIPROCITY  RECIPROCITY 


FROM 


2. 1.5-3  Cross-Axes  Coupling 


Cross-axes  coupling  is  the  motion  induced  in  a test  article  ortho- 
gonal to  the  direction  of  input  force  or  motion.  Figure  28  overlays  cross- 
axes coupling  impedance  on  an  axially  aligned  impedance  measurement  for  the 
same  accelerometer  location.  Accounting  for  cross-axes  coupling  was  not 
required  for  response  prediction^  on  this  building.  The  frequency  magnitude 
of  the  external  blast  loads  were  very  small  in  the  frequency  range  where 
cross-coupling  was  a significant  proportion  of  the  regular  impedance  measure 
ments.  The  root  sum  square  for  all  four  walls  contributed  less  than  5%  to 
vertical  motion.  In  the  present  ongoing  missile  silo  tests,  cross  axes 
coupling  is  significant  and  response  predictions  will  have  to  include  this 
effect . 

2.1.5-^  Symmetry 

Where  a test  article  or  structure  possesses  geometric  symmetry,  or 
at  least  reasonable-appearing  symmetry,  testing  may  be  shortened  by  measure- 
ments on  one-half  or  one-quarter  of  the  structure,  as  the  case  may  be. 
Decisions  to  use  symmetric  measurements  should  be  made  by  impedance  surveys. 
Where  unsymmetric  elements  or  components  exist,  either  local  measurements, 
substitution,  or  a strategy  of  redistribution  of  the  measurement  matrix  can 
often  be  judiciously  used.  Figure  29  shows  inertance  transfer  measurements 
overlaid  for  symmetry  comparison  on  a Ballistic  Missile  Defense  System 
Impedance  Test  Unit  (hog  trough). 


2.1.6  SYSTEM  FUNCTIONS 


Since  1972  numerous  impedance  projects  have 
ascertain  system  functions  and  to  predict  response  mot 
variety  of  external  threats.  From  these  projects,  the 
been  selected  because  of  project  uniqueness,  phenomena 
niques  employed.  The  data  presented  are  consistent  in 
factored  for  declassification. 


been  conducted  to 
ions  therefrom  for  a 
following  topics  have 
observed,  and  tech- 
this  repKjrt  but 
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2. 1.6.1  SAFEGUARD,  Perimeter  Acquisition  Radar  Building  (PARB) 


Impedance  measurements  were  made  upon  a I/12-scale  model  PARB  and 
the  prototype  PARB.  The  model  PARB  (Fig.  30)  is  located  at  the  Defence 
Research  Establishment,  Suffield,  Canada,  where  it  had  previously  been  sub- 
jected to  air-blast  and  ground-shock  loading  from  a 500-ton  TNT  explosion  in 
the  EVENT  DIAL  PACK.  The  prototype  PARB  (Fig.  31)  is  an  operational  radar 
system,  formerly  part  of  the  SAFEGUARD  system  located  in  North  Dakota.  The 
prototype  building  is  a reinforced-concrete  protective  structure  that  measures 
125  ft  in  height,  19k  ft  on  its  front,  and  210  ft  on  its  side,  with  exterior 
walls  3 ft  thick.  Pictures  of  the  vibrator  used  and  the  impedance  processing 
equipment  are  covered  in  Figures  32  and  33* 

Typical  systems  functions  of  inertance  magnitude,  phase,  and 
impulse  functions  are  given  in  Figure  3k  for  the  scale  model  and  in  Figure  35 
for  the  prototype.  These  nonparametri c functions  have  been  digitally  proc- 
essed into  final  form  for  use  in  response  predictions  of  the  structures. 

Response  predictions  for  a few  locations  involve  a large  number  of 
measurements  taken  over  the  structure  external  surfaces.  For  efficiency,  these 
predictions  can  be  used  to  predict  responses  at  additional  internal  locations 
of  the  building.  This  is  facilitated  by  use  of  transfer  functions.  The 
response  prediction  at  the  first  location  is  multiplied  by  a transfer  function 
to  another  location  to  obtain  the  response  at  the  new  location.  Transfer 
functions  are  obtained  by  the  complex  ratio  of  the  motion  at  a remote  location 
to  the  motion  at  the  drive  point  location,  as  illustrated  in  Figure  36. 

Transfer  functions  are  represented  as: 


t(u}) 


Calculation  of  responses  at  other  locations  in  the  structure  were  made  in 
accordance  with  the  following  equations: 


iij  (u))  “ Xg((D)  T.^j(a)) 


5k 


(Text  continued  on  p.  65.) 
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FIGURE  31.  PROTOTYPE  PERIMETER  ACQUISITION  RADAR  BUILDING  (PAKB) , NORTH  DAKOTA 
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FIGURE  32.  ELECTROHYDRAULIC  SHAKER  USED  FOR  MEASUREMENTS  OF  THE 
PROTOTYPE  PARS  (Waterways  Experiment  Station) 
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FIGURE  33.  INSIDE  VIEW  OF  INSTRUMENTATION  AND  IMPEDANCE  PROCESSING  VAN 
(Waterways  Experiment  Station) 
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FIGURE  35.  TYPICAL  PROTOTYPE  DIGITIZED  MEASUREMENT  IN  FINAL  FORM  FOR  COMPUTATIONS 


61 


FIGURE  36.  METHOD  OF  MEASUREMENT  OF  TRANSFER  FUNCTIONS 


where 


Frequency  response  at  location  i 

Frequency  response  calculated  for  center  of  the  5th  floor 
Transfer  function  location  i to  5th  floor 


X.(w) 

Xjlw) 


x.(t)  » 3-'|x.((u)| 

I III 

Measured  transfer  functions  to  the  Ath,  3rd,  and  2nd  floors  of  both 
the  nxDdel  and  prototype  structures  are  presented  in  Figures  37  and  38. 


Geometric  scaling  of  the  1/12-scale  model  to  the  prototype  is  given 
in  Table  1 . 


TABLE  1.  GEOMETRIC  SCALING 


Funct ion 

Prototype 

Mode  1 

Comp  1 i ance 

(^) 

' P 

Mob i 1 i ty 

1 nertance 

(f) 

P 

1 /x  \ 

T72?  \f/ 

' m 

F requency 

TT 
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Scaling  the  model  data  geometrically  did  not  provide  a good  fit. 
Several  not-to-scale  effects  existed  in  the  prototype  that  had  not  been  or 
could  not  be  accounted  for  In  the  model.  These  effects  included: 

a.  Structure 

• Prototype  roof--radar  base  slab 

• Prototype--2nd  floor  mezzanine 

• Prototype--ut i I i ty  tunnel 

• Model --excess i ve I y large  base  slab 

b.  Equipment  in  Prototype  Structure 

• Tactical  support  equipment  on  shock- i so  1 a ted  platforms, 

2nd  floor  365,000-lb  sprung  weight  of  22  platforms,  and 

3rd  floor  203,000-lb  sprung  weight  of  k platforms 

• SAECA  communication  equipment  on  shock- i sol ated  platforms, 
3rd  floor  86,126-lb  sprung  weight  of  3 platforms 

• Weapons  systems  equipment  hard  mounted  on  3rd  floor 

The  radar  base  slab  on  the  prototype  roof  stiffens  the  area  covered 
and  adds  weight.  The  2nd  floor  mezzanine,  the  utility  tunne'  of  the  prototype, 
and  the  extra  large  base  slab  in  the  model  contribute  to  the  scaling  problem. 
Equipment  installed  in  the  prototype,  particularly  on  the  2nd  and  3rd  floors, 
has  a decided  effect  in  altering  the  floor  frequency  responses.  The  weapon 
system  equipment  hard  mounted  to  the  3rd  floor  is  massive  and,  in  the  case  of 

the  precise  power  units,  acts  as  stiffeners.  These  power  units  cover  approxi- 

mately one-third  of  the  3rd  floor  area.  The  numerous  shock-isolated  platforms 
are  dynamically  responsive  over  a wide  frequency  band  and  dynamically  inter- 
act with  the  floors  to  alter  floor  responses  substantially. 

t A scaling  comparison  in  both  magnitude  and  frequency  was  made 

between  representative  and  paired  inertance  measurements  of  the  model  and 
prototype  buildings.  These  functions,  such  as  those  shown  in  Figures  3^  and 
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35,  are  compared  for  the  5th  floor  and  similarly  for  each  floor.  Scaling 
parameters  were  determined  by  keying  corresponding  model  and  prototype  func- 
tions at  their  peak  amplitudes  and  then  iteratively  determining  amplitude 
and  frequency  scale  factors  for  a minimum  error  difference.  The  scale  factors 
determined  by  this  procedure  are  given  in  Table  2 for  each  location. 

TABLE  2.  SCALING  FACTORS.  MODEL  TO  PROTOTYPE,  FOR  GEOMETRIC 
AND  INERTANCE  MEASUREMENT  SCALING 


Scale 

Factors 

Oi f fercnce 

Inertance 

Geometric 

(Reference  Geometric) 

Floor 

Oi splacement 
Factor,  a 

Frequency 
Factor,  b 

Di spl acement 
Factor,  A 

F requency 
Factor,  A 

Displacement , 

% 

Frequency , 

5th 

18.1. 

12 

12 

♦53 

-15 

l*th 

10.5 

12 

12 

-12 

+6 

Jrd 

25.3 

12 

12 

+ 1 I 1 

+ 17 

2nd 

3A.5 

■1 

12 

12 

*188 

-19 

AA  7696 


2. I. 6. 2 System  Functions  of  Massive  Equipment 

System  functions  were  measured  on  ten  large  shock- i sol ated  platforms, 
which  ranged  in  weight  from  7,500  lb  to  21*4,000  lb.  The  method  of  measurement 
presented  herein  is  of  more  interest  than  the  systems  measured,  for  its  general 
applicability  to  in-place  massive  machinery  and  equipments.  Typical  shock- 
isolated  platforms  for  which  transfer  functions  were  measured  are  pictured 
in  Figure  39- 


The  indirect  method  of  transfer  function  measurement  is  illustrated 
in  Figure  40,  which  measures  the  transfer  function  at  a point  of  practical 
and  phy.ical  convenience.  The  size  of  machines  and  equipment  often  make  it  • 

impossible  to  jack  up  an  equipment  mounting  leg  and  insert  a vibrator.  By 
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(a) 


FIGURE  39.  TYPICAL  SHOCK  ISOLATED  PLATFORMS  UPON  WHICH  TRANSFER  FUNCTIONS 
WERE  MEASURED 
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(b) 


FIGURE  39.  (CONCLUDED) 
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FIGURE  40.  PLATFORM  PARPP-D:  INDIRECT  TRANSFER  FUNCTION  MEASUREMENTS 
t(u)  - F.U)/F,(u))  • (u)) 


measuring  in  a direction  opposite  the  input  — from  the  equipment  to  its  mount- 
ing points--a  force-output  to  force-input  ratio  is  taken.  This  ratio  is 
identical  to  the  acceleration  ratio  required  for  the  input  motion  direction. 
From  the  vibration-generator  position  to  other  locations  on  the  equifjment, 
acceleration  ratios  are  used.  The  indirect  transfer  function  is  composed  of 
two  parts,  a complex  ratio  of  output  force  to  input  force  and  a complex  ratio 
of  output  acceleration  to  input  acceleration.  These  two  ratios  are  multiplied 
together  to  form  the  transfer  function  from  each  isolator  input  to  a specific 
platform  location.  The  computations,  which  are  required  for  each  drive  point, 
are  expressed  as  follows: 

Force  ratio  per  drive  point 
|F 

' O ] 

where 

F (u))  Complex  output  force  at  equipment/building  junction 

o 

Pj  (uj)  ■ Complex  input  force  at  drive  point 
9,4i  ■ Phase  angles 

Acceleration  ratio 


Tp(li))  . 


F^(u.) 
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where 

Xo(a.) 
*j  (u)) 

6.Y 


Complex  platform  accelerometer  signal 
Complex  drive-point  accelerometer  signal 
Phase  angles 


Transfer  function 


t(u>)  - Tp(u))  • T^(u)) 


*0 


(u)) 

U) 


The  transfer  function,  computed  as  a product  of  both  the  force  and 
acceleration  ratios,  is  summed  for  all  isolators  of  a platform  to  generate  an 
overall  (global)  platform  transfer  function. 

A plot  of  transfer  function  magnitude  (Fig.  ^1)  reveals  that 
resonance  amplifications  from  modes  occur  from  0.5  Hz  to  10  Hz,  with  magni- 
tudes from  1.2  to  2.3-  The  resonance  at  17  Hz  for  a magnitude  of  0.7  was 
identified  during  measurements  as  spring  surge  in  the  coil  isolators. 

Another  example  of  a transfer  function  is  shown  (Fig.  k2)  for 
Platform  D.  This  platform  measures  A ft  x 12  ft,  weighs  13.100  lb,  and  is 
supported  by  four  undamped  coil-spring  isolators,  pendulum  mounted.  Three 
a.c.  switchboard  cabinets  are  mounted  on  this  platform.  Data  measurements 
are  from  35  Hz  to  500  Hz,  and  rigid-body  modes  have  been  added  to  Figure  k2 
for  vertical,  pitch,  and  roll  frequencies  at  0.85,  1.08,  and  1.^5  Hz. 

A comparison  is  made  in  the  plot  of  Figure  42  between  the  measured 
resonance  frequencies  and  the  modal  frequencies  of  a finite  element  model  for 
this  platform.  It  is  to  be  noted  that  28  modal  frequencies  were  found  by  the 
model  and  approximately  77  modal  frequencies  are  present  by  measurement. 


\ 
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FIGURE  1*1.  PARPP  CONT 


2. 1.6.3  Buried  Structures 


A series  of  experiments  performed  by  the  Waterways  Experiment  Station 
{Ref.  37)  in  Project  ESSEX  was  made  to  assess  the  effects  of  soil  coupling  on 
structures.  Two  reinforced  box  structures  were  constructed  with  both  struc- 
tures having  the  same  interior  dimensions  of  16  ft  x A ft  x ^ ft.  The  struc- 
tures differed  in  wall  thicknesses  of  13  in.  and  5-6  in.  The  results  of 
impedance  tests  on  the  buried  structures  indicated  (Figs.  43  and  44)  that  the 
effect  of  soil  cover  on  the  resonant  frequencies  and  damping  appear  to  be 
dependent  upon  the  structure/soil  impedance  ratio.  Further  work  is  planned 
by  which  the  soil  effect  will  be  extracted  and  modeled  from  the  experimental 
data.  Extraction  procedures  would  be  in  accordance  with  the  method  presented 
on  component  interconnections  in  Reference  9- 

2. 1.6.4  Electrical  Distribution  Center,  SAFEGUARD  System 

The  Electrical  Distribution  Center  of  the  SAFEGUARD  SYSTEM  is  a 
reinforced-concrete  protective  structure  flush  to  ground  surface  (Fig.  4$). 

This  structure  measures  26  ft  x 43  ft  and  is  IS  ft  deep.  The  ceiling  is 
approximately  2 ft  in  thickness.  Internal  electrical  equipment  is  mounted 
on  a shock- i sol ated  platform  (30  ft  x 20  ft)  that  weighs  on  the  order  of 
40,000  lb.  The  platform  is  supported  at  each  corner  by  pneumatic  isolators, 
pendulum  mounted  to  the  ceiling  (Fig.  46). 

Measurements  made  at  the  attachment  location  of  the  shock  isolator 
with  and  without  isolator  attached  (Fig.  47)  show  differences.  This  measure- 
ment difference,  although  small,  illustrates  the  sensitivity  of  the  measure- 
ment system  even  for  components  having  a decidedly  large  impedance  mismatch. 
Extraction  of  the  isolator  input  impedance  will  be  made  during  the  data 
analysis  phase  of  this  project.  The  extraction  method  will  be  in  accordance 
with  procedures  presented  on  interconnection  in  Reference  35.  Figure  27 
also  shows  nonlinear  reciprocity  measurements  on  this  structure  between  roof 
and  floor. 


(Text  continued  on  p.  80.) 
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TYPICAL  IMPEDANCE  FUNCTIONS  FOR  THICK  WALLED  (13  IN.)  BOX  STRUCTURE 
(16  FT  X l(  FT  X A FT)  SHOWING  EFFECTS  OF  SOIL  COVER 
(Courtesy  Waterways  Experiment  Station) 
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FIGURE  kk. 


TYPICAL  IMPEDANCE  FUNCTIONS  FOR  THIN  WALLED  (S 
(16  FT  X A FT  X A FT  INTERIOR)  SHOWING  EFFECTS 
(Courtesy  Waterways  Experiment  Station) 


6 IN.)  BOX  STRUCTURE 
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FIGURE  A5. 


ELECTRICAL 

PROTECTIVE 


DISTRIBUTION  CENTER,  A FLUSH-TO-GROUNO  SURFACE  REINFORCED 
STRUCTURE  (Two  electrohydraul ic  vibrators  on  roof) 
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FIGURE  kj.  DRIVING  POINT  IMPEDANCE  TESTS  OF  ELECTRICAL  DISTRIBUTION  CENTER 
AT  ATTACHMENT  LOCATION  OF  A SHOCK  ISOLATOR  (Measurements  made 
with  and  without  isolator  attached) 
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2. 1.6. 5 Other  Projects 


BMP  Test  Structure 

Impedance  measurements  shown  in  Figure  3 for  chirp  tests  and 
Figure  12  for  random  tests  were  made  on  the  Ballistic  Missile  Division 
(Bmp)  test  structure  pictured  in  Figure  48.  This  structure  of  reinforced 
concrete  was  used  as  a test  device  to  study  various  methods  of  vibration 
excitation  and  data  processing  to  improve  accuracy  of  data  and  to  reduce  test 
and  data  processing  costs. 


Sprint  Missile  and  Silo 

Impedance  measurements  of  Figure  23  exhibit  effects  of  nonlinear 
soil  coupling  to  a missile  silo  wall.  These  measurements  were  recently  taken 
on  an  ongoing  impedance  project  for  the  Sprint  Silo  and  Sprint  Missile 
(Fig.  49). 


Hangeret  te 

Hangerettes  (Fig.  50)  are  used  to  protect  aircraft  from  both  con- 
ventional bombs  and  nuclear  threats.  Mobility  measurements  (Vel/Force)  were 
made  (Fig.  51)  from  the  inside  surface  of  the  re i nforced-concre te  arch  to  the 
external  soil  surface.  These  measurements  in  a reciprocity  sense  assume 
linearity  from  which  response  predictions  can  be  made,  as  will  be  covered  in 
Section  2.2. 


(Text  continued  on  p.  85.) 
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FIGURE  48.  BALLISTIC  MISSILE  DIVISION  (BMD)  IMPEDANCE  TEST  STRUCTURE 
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FIGURE  50.  HANGARETTE-AI RCRAFT  SHELTER  SHOWING  LOCATION  OF 
PRESSURE  AND  VELOCITY  GAGES  MIXED  COMPANY  EVENT 
(500-T0N  TNT  TEST) 
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(b)  Mobility  impulse  function 


FIGURE  51.  MOBILITY  OF  HANGARETTE  FROM  MID-ARCH  TO  SOIL  SURFACE  (Courtesy 
Waterways  E;<periment  Station) 
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2.2  PREDICTION  OF  SYSTEM  RESPONSES  USING  MEASURED  IMPEDANCE  AND 
TRANSFER  FUNCTIONS 


2.2.1  RESPONSE  OF  ABOVEGROUND  STRUCTURES  TO  AIR  BLAST 

The  method  used  to  calculate  the  response  accelerations  of  the 
1/12-scale  model  and  prototype  Perimeter  Acquisition  Radar  Buildings  (PARB) 
is  similar  for  each.  Verification  or  checks  on  the  accuracy  of  the  response 
predictions  is  made  by  comparing  the  predictions  with  response  records  of  the 
scale  model  taken  during  Event  DIAL  PACK  (Ref.  3).  Event  DIAL  PACK  was  a 
500-ton  TNT  explosion  that  subjected  air-blast  and  ground-shock  loads  to 
several  targets.  alculations  for  both  structures  were  made  for  the  vertical 
direction  only,  and  as  such,  the  traveling  air-blast  pressure  loading  over 
the  roof  surface  comprises  the  external  forcing  function.  Exclusion  of  the 
cross-coupling  contributions  to  structural  motions  from  pressure  loads  on 
the  external  four  walls  is  considered  to  have  minimal  effect  on  the  calculated 
responses.  As  discussed  in  Section  2.1,  the  rapid  rolloff  of  pressure  ampli- 
tudes with  frequency  permitted  this  reduction  in  computational  effort. 

All  data  and  calculations  are  internally  consistent  and  show  com- 
parable amplitudes  for  pressures,  motions,  yield,  time,  and  frequency. 

However,  the  absolute  values  have  been  normalized  so  that  this  report  could 
be  declassified. 

2. 2. 1.1  1/1 2-Scale  Model 

Response  predictions  use  inertance  functions,  a pressure  zone  for 
each  inertance  function,  air-blast  pressure  functions,  pressure  area  zone 
engulfment,  traveling  pressure  wave  function,  and  transfer  functions.  Each 
of  the  above  elements  is  described  and  quantified  in  the  following  paragraphs. 
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2.2.  1.1.1  Tncrtan^^  Functions  and  Associated  I^essure  Area  Zoning 

The  use  of  inertance  measurements  to  predict  the  response  of  a 
structural  system  to  blast  loading  is  basically  a special  application  of 
the  impedance  technique.  An  inertance  test  measures  the  point  acceleration 
response  of  a structure  due  to  a point  force  load.  Blasts  such  as  the 
air-pressure  loads  acting  on  the  model  PARB  structure  in  Event  DIAL  PACK  is 
a continuous  forcing  function  acting  on  a continuous  structure.  Thus,  the 
sampled  Impedance  information  must  be  summed  over  the  structure  (as  an 
approximation  of  an  integral)  to  obtain  the  internal  point  responses  of  the 
structure  due  to  loads  acting  everywhere  on  the  external  surfaces  of  the 
structure. 


Each  inertance  function  must  be  associated  with  a specific  area 
of  the  structure  upon  which  the  air-blast  pressure  acts.  This  assumes  that 
the  inertance  function  is  essentially  constant  for  this  area.  Finite-element 
plate  models  were  initially  used  to  map  the  model  and  to  provide,  thereby, 
preliminary  measurement  and  drive  point  locations.  Additional  refinement 
of  drive  point  locations  was  made  during  field  testing  where  quick-look 
analog  impedance  measurements  verified  that  a reasonable  density  of  measure- 
ments for  the  model  structure  would  be  acquired.  The  smooth  change  in 
functions  obtained  when  moving  from  one  area  to  the  next  indicated  that  a 
reasonable  number  of  measurements  consistent  with  computational  costs  had 
been  obtained.  Potentially,  interpolation  of  the  measurements  could  have 
been  made  during  computation,  should  the  need  have  arisen.  Additional  measure 
ments  also  showed  a high  degree  of  symmetry  for  the  structure  (see  Sec.  1), 
and  this  fact  permitted  a higher  density  of  measurement  per  surface  area. 

Boundaries  for  local  areas  of  each  acceptance  measurement  were 
drawn  by  observation  and  evaluation  of  the  data.  For  the  29  measurements 
selected,  the  associated  impedance  zone  and  input  pressure  zone  for  each 
measurement  is  shown  in  Figure  52. 
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FIGURE  52.  ZONED  ROOF  ON  PARB  MODEL  SHOWING  IMPEDANCE  ZONES,  AIR-BLAST  ZONES, 
AND  IMPEDANCE  MEASUREMENT  DRIVE  LOCATIONS 
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2 . 2 . 1 . 1 . 2 Prcesurt’  Input  Developrrx’nt 


Only  three  air-blast  measurements  (Ref.  29)  on  the  roof  were 
obtained  from  Event  DIAL  PACK  (Figs.  53  <snd  5^)-  The  aboveground  structire 
perturbs  the  air-blast  loading  by  overexpanding  the  transient  air  flow,  as 
is  particularly  observable  in  the  pressure  notch  in  the  time  history  of 
Figure  53b  and  in  the  frequency  spectrum  of  Figure  5^b.  Hence,  test  data 
rather  than  conventional  blast  curves  (Figs.  53d  and  5^d)  were  used  as 
forcing  functions. 

The  air-blast  zones  shown  in  Figure  52  were  configured  on  a simple 
geometric  basis.  The  air-blast  functions  were  considered  to  be  constant 
within  each  air-blast  zone.  The  potential  for  extrapolation  and  interpolation 
exists  for  this  type  of  data. 

a.  Pressure  Area  Zones 

For  a local  surface  area,  A.|^,  which  is  associated  with  one  transfer 
inertance  function,  it  is  necessary  to  account  for  the  velocity  of  the  air 
load  as  it  traverses  the  surface.  This  accounting,  as  illustrated  in 
Figure  55,  is  required  to  assure  a proper  surface  loading. 

Upon  arrival  of  the  air  blast  at  time  x.,  the  air  blast  successively 
loads  the  surface  in  increments  AX  where  nAX  = d.  the  length.  The  incre- 
mental loading  time  is  given  by:  At  = AX/V.  This  incremental  loading  subdi- 
vides the  local  surface  into  segments  A..  (AX/d.).  The  fineness  of  AX  is 

I k I 

obtained  from  the  At  selected  and  can  be  as  small  as  the  digitizing  rate. 
Depiction  of  air-blast  engulfment  and  the  development  of  an  effective  pressure 
as  seen  by  a local  area  is  shown  in  Figure  56. 


(Text  continued  on  p.  93.) 
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FIGURE  53-  AIR-BLAST  MEASUREMENTS  ON  MODEL  PARB  ROOF  AND  CONVENTIONAL  AIR-BLAST 
REPRESENTATION  OF  PRESSURE-TIME  HISTORIES 
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OF  FOURIER  MAGNITUDE  PRESSURE  SPECTRA 
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FIGURE  55-  MOVING  AIR  LOAD  OVER  A LOCAL  SURFACE  AREA  ASSOCIATtt  . * 
TRANSFER  INERTANCE  FUNCTION  TO  AN  INTERNAL  POINT  IN 
STRUCTURE 
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A pressure-t ime  history  of  frontal  velocity  V, 
translated  and  sectioned  over  the  surface  of  length 
increments  is  given  by: 


which  is  successively 
d.  equal  to  nAX 


Peff(t) 


I \ ^ I.  I,  ^ I.  2AX>  , 2AX> 

p(t)  + p(t  - u(t  - * p(t  - -rr-l  u(t  - — ) 


...  p(t  - 

and  transformation  in  the  frequency  domain  gives: 


(n  - 1)AX  , (n  - 1)AX 
V ■ V 


(1) 


p(.) 


1 + c 


.AX 

-j— 


+ e 


.2AX 
•J—  ^ 


.3AX  .(n-l)AX 

“ *J w “I 

+ e ■ ...  + ''  1 (2) 


and  the  expression  of  Equation  2 converges  to: 


P^ff(w) 


P(u.) 


. /nAX\ 
sin 


n sin 


,AX. 


. (n-1) AX 
J 2V 


(3) 


for  small  angles  (6  - sin  9)  the  engulfing  function  of  Equation  3 may  be 

represented  by: 


s I n 


,nAX> 


sin  (27)“ 


sin 


/ nAX  \ 


,AX. 

n 


(M 


which  is  the  Sj  (X)  function  in  the  frequency  domain  such  as  represented 
in  Figure  57a  and  an  approximate  box  car  function  in  the  time  domain.  The 
cyclic  nature  of  the  actual  engulfing  function  (Eq.  3)  is  illustrated  in 
Figure  57b  and  implies  a critical  selection  of  a time  step  At  (dependent 
on  mesh  length/wave  velocity)  for  the  frequency  limit  involved.  Somewhat  akin 
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(•)  Cn9ulfln9  function  Mitn  ci««  step  idocCM  for  proper  frMupnev  linic  r«Quirtfn«nc» 


(b)  Engulfing  function  with  time  ttpp  improperly  teleetep  for  frequency  limit  requi rmnents 


FIGURE  57.  PLOT  OF  MAGNITUDE  OF  ENGULFING  FUNCTION  FOR  TWO  DIFFERENT 
TIME  STEPS  At  - (MESH  LENGTH/BLAST  VELOCITY) 


94 


to  the  Nyquist  folding  frequency,  the  engulfing  function,  as  shown  in 
Figure  57a  has  a tinoe  step  At  (dependent  on  mesh  length/velocity)  and  a 
frequency  cutoff  1/2  At  (0.1  msec)  for  5 KHZ.  These  values  are  more  than 
adequate  for  the  model  inertance  function  where  measurements  were  made  to 
3 KHz.  Obviously,  serious  effects  are  encountered  when  too  large  a time  step 
(mesh  length/velocity)  is  taken  in  using  a traveling  wave,  since  the  results 
will  ultimately  include  the  effects  of  the  function  as  shown  in  Figure  57b, 
where  frequency  amplitudes  above  cutoff  are  contained  in  the  data  or  computa- 
tions. Not  too  obvious  in  the  engulfing  functions  is  the  maximum  amplitude 
in  the  frequency  function  of  unity  and  the  area  of  the  approximate  boxcar 
function  in  the  time  domain  of  unity.  Figure  58  provides,  for  a traveling 
load  over  a surface,  the  maximum  mesh  length  at  the  folding  frequency  for  a 
range  of  traversing  velocities. 


b.  Traveling  Wave  over  a Surface 

Somewnat  similar  to  the  engulfing  function  for  local  areas  (mesh 
size)  is  the  wave  travel  over  the  entire  structure.  In  this  case,  the  travel 
distance  or  time  increment  is  from  the  boundary  of  one  local  area  to  the 
next.  Where  the  boundary  distances  are  equal  and  the  blast  velocity  constant, 
the  traveling  wave  function  would  be: 


where 


. /NAt> 
sin  (-^)u) 

JaT 

sin  (-j-)u) 


. (N-l)At 
■'  2 


N • Number  of  equal  local  areas  on  roof 
At  ■ Traverse  time  across  local  area 


(5) 
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As  before  for  small  angles,  d 
becomes : 


sin  3 and  the  function  of  Equation  5 


sin  (-y)u 
/ it  , 


■ (N-i)i:it 

J 2 


(6) 


The  traveling  wave  for  the  roof  is  a summation  having  a maximum 
amplitude  in  the  frequency  domain  of  NAt  (the  total  traverse  time  of  the 
wave  over  the  roof).  The  amplitude  of  the  approximate  boxcar  function  is 
unity  and  has  an  area  of  NAt  in  the  time  domain. 


The  local  areas  of  the  roof  are  not  equal,  and  the  foregoing  dis- 
cussion is  useful  only  in  understanding  the  characteristics  of  the  function. 
The  unequalness  of  the  local  roof  areas  requires  that  the  function  of 
Equation  5 be  placed  in  the  form 


-jt  Ui  -jTjUl 

)+e  t-e  +e  ■’+...  e 


(7) 


where  t.  • time  of  arrival  at  each  local  area  boundary 
1 

2. 2. 1.1. 3 Beeponae  Prediction  at  First  Location 

Response  calculations  were  made  in  the  vertical  direction  from  roof 
loads  of  a traveling  air  blast  for  the  center  of  the  5th  floor  of  the  1/12-scale 
model.  The  method  of  calculation  Is  given  as: 


Response 


Local 
Inertance 


Loca  I 
Area 


Local 

Pressure 


Local  \ /Traveling 
x Engulfmentj  I Wave 
Function  / 'Function 


Computations  were  performed  in  the  frequency  domain  as  represented  by 
Equation  8 and  transformed  to  the  time  domain  by  Equation  9. 
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i th  inertance  function  [X/F.(u))] 
i th  local  area  (mesh)  of  external  load 
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i th  local  engulfment  function 
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■ Traveling  wave  function 
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Arrival  time  of  blast  wave  at  local  area  boundary 
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2.  2.1.1.  A Reeponae  Prediction  to  Other  Locations 

Response  predictions  at  additional  internal  locations  (i.e.,  loca- 
tions in  addition  to  the  "first")  of  the  building  are  facilitated  by  use  of 
transfer  functions.  The  response  prediction  at  the  first  location  on  the 
5th  floor  is  multiplied  by  a transfer  function  to  another  location  to  obtain 
the  response  at  the  new  location,  as  discussed  In  Section  2.1. 
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2.2. 1.2  Prototype  PARS 

Response  predictions  of  the  prototype  PARB  follow  the  same  procedure 
as  that  used  for  the  model.  Each  variable  will  be  discussed,  with  emphasis 
placed  where  differences  in  procedure  occur. 

Selection  of  local  impedance  zones  and  pressure  areas  for  each 
inertance  function  followed  the  methods  employed  for  the  model.  The  zones 
for  the  36  impedance  measurements  selected  are  shown  in  Figure  59-  Impedance 
measurements  A-1  and  A-29  were  used  on  the  basis  of  symmetry  in  two  other 
locations  (Fig.  59),  since  poor  data  records  were  obtained  in  Zones  3^  and  36. 

2.2. 1.2. 1 Preesure  Input  Development 

Pressure-time  histories  from  Event  DIAL  PACK  were  scaled  geomet- 
rically for  computations  on  the  prototype  PARB.  Pressures  for  both  model 
and  prototype  are  identical  but  scaling  of  characteristic  time  is  required. 


P (u>)  - X P (Xoi)  (13) 

p m 

Effects  on  waveforms  of  the  DIAL  PACK  pressure  records  by  scaling 
(Reynolds  number.  oVL/u)  to  the  prototype  were  evaluated  to  determine  whether 
serious  errors  would  result.  Information  on  transient  blast  waves  interacting 
with  structures  is  quite  limited.  Some  information  having  applicability  to 
the  PARB  configuration  has  been  found  and  evaluated  for  scale  factor  effects. 
These  data  are  from  shock  tubes,  water  tables,  and  nose  cones. 
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FIGURE  59-  ZONED  ROOF  OF  PARB  PROTOTYPE  SHOWING  IMPEDANCE  ZONES,  AIR-BLAST 
ZONES,  AND  IMPEDANCE  MEASUREMENT  DRIVE  POINTS  (DIRECT  AND 
RECIPROCITY) 
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2. 2. I. 2. 2  Aiv-Blast-Uave  Engulfment 


a.  Pressure  Area  Zones.  The  engulfing  function  developed  for  the 
model  (Eq.  was  applied  also  for  the  local  pressure  areas 
(impedance  zones)  of  the  prototype.  The  only  change  required 
was  a coarser  mesh,  owing  to  the  lower  cutoff  frequency  of 

300  Hz  (Fig.  58) . 

b.  Traveling  Wave  over  a Surface.  The  traveling  wave  function, 
Equation  7,  used  for  the  model  was  applied  to  the  prototype. 

2. 2. 1.2. 3 Eesponae  Prediation  at  First  Location 

Response  predictions  for  the  prototype  were  made  with  Equations  8 
and  9 in  a manner  similar  to  that  used  for  the  model. 

2. 2. 1.2.1*  Response  Predictions  for  Other  Locations 

Response  predictions  for  the  prototype  at  other  locations  were  made 
in  a manner  similar  to  that  used  for  the  model. 

2.2. 1.3  Development  of  Scaling  Relationships 

Scaling  relationships  for  the  1/12-scale  model  calculations  and 
DIAL  PACK  event  records  of  the  nwdel  PARS  are  developed  in  this  section. 

These  scaling  relationships  provide  for  both  geometric  scaling  and  scaling 
based  upon  comparison  of  inertance  functions  from  the  model  and  the  prototype 
from  Table  2. 

2. 2. 1.3.1  Geometric  Scaling 

a.  Areas 

Surface  area  of  the  air  blast  for  the  model  scaled  to  the  prototype  is 

A - A (I**) 

p m 


100 


where 


A “ Area  (sq  ft) 
p ” Prototype 
m « Model 

A » Scale  factor  (12) 


b.  Air  Blast 

The  air-blast  pressures  from  DIAL  PACK  are  scaled  for  the  prototype 
where  pressure  amplitudes  for  model  and  prototype  are  equivalent,  but  the 
time  duration  is  scaled  by  A(12). 


Pp(t) 


Pp(u)  • X Pi^(Xu)) 


(15) 


c.  Inertance  Functions 

For  simplicity,  a generalized  single-degree-of-freedom  model  was 
used  to  develop  the  scaling  relationship 


(16) 

(17) 

where 

F ■ Force  (lb) 

Mp  ■ Mass  (lb-sec  ^/in.)  - x\ 

K ■ Stiffness  (Ib/in.)  ■ XK 

p m 

2 

C ■ Damping  (ib-sec/in.)  - X C 

P m 


I"  - i)  -i! 


tV  (u) 


(U)  M - K)  -juiC 


\ « Scale  factor  (12) 

..  2 
X ••  Accelerations  (in. /sec  ) 

.r/  (u)  » Inertance  function  (in. /sec  -lb) 

Substitution  in  Equation  17  gives: 

,r/  (u)  ” ,<J  (Xiij)  (l8) 

P ^3  m 


d.  Geometrically  Scaled  Response 


Xp(w) 

• 2 r ' Si.(X(ij)  e ' ( 

i-o(x^  ' 1 ' ' '1  I 

(19) 

■’m 

X {(j) 
p 

« X (A(j) 

m 

' (20) 

X (t) 

P 

- - X m 

X m\X/ 

(21) 

where 


P 


m 

X 

tX/j  (Xai) 
A. 

I 

P.  (Xiu) 


S i . (Xu)) 


- j T . X(i) 
e ' 


Prototype 
Mode  1 


Scale  factor 

i th  impedance  function  [X/F.(u)] 
i th  local  area  of  external  load 
ith  static  external  pressure  load 


i th  local  engulfment  function 


sin  (n^) 
n sin 


) 


Traveling  wave  function 


+ e 


•jitXii)  ^ -2j:ltu; 


Arrival  time  of  blast  wave  at  local  area  boundary 
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2.2.1. 3- 2 Scaling  by  Use  of  Measured  Inertance  Functions 


a.  Scaling  of  Areas,  Pressures,  and  Wave  Engulfment 

Scaling  of  areas,  pressure,  and  wave  engulfment  are  geometric  as 
presented  in  Section  2. 2. 1.3.1. 


b.  Inertance  Function  Scaling 

Several  not-to-scale  conditions  existed  between  the  prototype  and 
model  PARS.  These  conditions  cannot  really  be  quantitatively  established, 
although  discussion  and  identification  of  each  is  given  in  Section  1.  A 
scaling  comparison  in  both  magnitude  and  frequency  was  made  between  rtpre 
sentative  and  paired  acceptance  measurements  of  the  model  and  prototype 
buildings.  Scaling  parameters  were  determined  by  keying  corresponding  model 
and  prototype  functions  at  their  peak  amplitudes  and  then  iteratively  deter- 
mining amplitude  and  frequency  scale  factors  for  a minimum  error  difference. 
The  scale  factors  determined  by  this  procedure  are  presented  in  Tables  2 and  3 
for  each  locat i on . 

c.  rteasured  Inertance  Scaling  (not-to-scale  effects) 


N ( 

i 1 

X (a.) 

" S 

y Jl.  (buj) 

X^A.  • XP.  (Xu) 

y 

i-ol 

!ab 

I 1 

•Xp(.) 

x3 

,2 

X (bw) 
m 

ab 

X (t) 

. _xi 

X (f) 

P 

ab^ 

m'  b 1 

Si  . ( 'uj) 

I 


(22) 

(23) 

(2k) 
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TABLE  3.  ACCELERATION  RESPONSE  SCALING,  MODEL  TO  PROTOTYPE,  OF  GEOMETRIC 
AND  INERTANCE  MEASUREMENT  SCALE  FACTORS 
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2. 2.1. A Comparison  of  Prototype  Responses  to  Response  Scaled  from  Model 

Acceleration  response  motions  from  DIAL  PACK  records  and  model 
calculations  are  scaled  to  the  prototype  calculations  and  presented  in 
Figures  60  through  67.  Scaling  was  made  by  inertance  function  ratios,  since 
geometric  scaling  did  not  provide  a good  fit. 

Overall  similarity  of  scaled  DIAL  PACK  and  model  calculations  to 
prototype  calculations  is  guite  good  for  acceleration-time  histories  and 
Fourier  magnitude  freguency  spectra. 

Deviations  noted  in  the  calculated  rise  and  fall  times  for  time 
histories  involve  both  data  measurements  and  subsequent  processing.  Some 
improvements  in  technique  were  applied  to  the  prototype  calculations.  Differ- 
ences observable  in  the  Fourier  fansforms,  particularly  in  the  Ath  and 
3rd  floors,  are  attributable  to  not-to-scale  effects  in  the  prototype  and  to 
measurement  processing  techniques  for  f loor-to-f loor  transfer  functions. 


(Text  continued  on  p.  113.) 


COMPARISON  OF  PROTOTYPE  AND  SCALED  MODEL  PARS  ACCELERATION-TIME 
HISTORIES:  ATH  FLOOR  CENTER 


HISTORIES;  3RO  FLOOR  CENTER 


I 


HISTORIES:  2ND  FLOOR  CENTER 


COMPARISON  OF  PROTOTYPE  AND  SCALED  MODEL  FOURIER  MAGNITUDES 
OF  ACCELERATION:  5TH  FLOOR  CENTER 


OF  ACCELERATION:  4TH  FLOOR  CENTER 


OF  ACCELERATION:  3RD  FLOOR  CENTER 


OF  ACCELERATION:  2N0  FLOOR  CENTER 


2.2.2  RESPONSE  OF  MASSIVE  EQUIPMENT  TO  SPECIFIED  BASE  MOTION  INPUT  CRITERIA 


2.2.2. 1 Environmental  Threat 

The  shock  environments  Internal  to  structures  at  the  mounting 
locations  of  equipments  were  specified  as  undamped  shock  spectra.  The 
environmental  shock  spectra  selected  to  test  and  predict  the  response  of  all 
weapons  systems  In  protective  structures  were  those  designated  as  Annex  M 
and  shown  in  Figure  68a.  These  worst-case  (envelope)  environments  were  chosen 
because  they  equal  or  exceed  the  envelope  of  local  input  shock  spectra  for  all 
structures  of  the  system  at  practically  all  frequencies,  and  are  therefore 
considered  the  reference  shock  spectra. 

For  analysis  and  test,  a time  history  (acceleration,  velocity,  or 
displacement,  versus  time)  is  required  as  input  to  predict  and  measure 
platform  motion.  However,  the  specified  envelope  spectrum  (Annex  M)  is 
generalized  to  envelop  numerous  threat  histories.  Such  lack  of  uniqueness 
therefore  requires  that  many  tests  or  calculations  be  made.  An  alternative 
was  developed  by  Yang  (Ref.  38)  to  generate  a synthesized  time  history  that 
populates  specified  shock  spectra  for  the  protective  structures  and  thereby 
assures  a conservative,  or  worst-case,  test.  This  synthesized  time  history 
is  the  summation  of  a series  of  oscillations  at  selected  frequencies  over 
the  range  of  the  shock  spectra.  Oscillations  at  each  frequency  are  from 
1-1/2  oscillations  (3  one-half  cycles)  to  6-1/2  oscillations  (13  one-half 
cycles).  The  amplitude  of  these  oscillations  is  constrained  to  a sine-beat 
envelope.  Synthesized  waveforms  that  populate  the  input  shock  spectra  of 
Annex  M are  shown  in  Figure  68b. 

2. 2. 2. 2 Response  Predictions 

The  motions  of  the  10  shock- i sol ated  platforms  were  calculated  with 
input  from  both  the  overall  Annex  M motion  and  the  specified  local  motion  for 
each  platform  location.  In  the  frequency  domain,  the  Input  functions  are 
multiplied  by  the  system  functions,  followed  by  an  Inverse  Fourier  transforma- 
tion to  obtain  the  platform  acceleration-time  history.  In  the  time  domain, 
input  time  history  of  the  building  is  convolved  with  the  platform  impulse 
function. 
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(b)  Synthesized  waveform,  vertical  acceleration 
time  history  to  match  Annex  M shock  spectrum 


FIGURE  68.  VERTICAL  OVERALL  (ENVELOPE)  BUILDING  INPUT  MOTION  FOR  ALL 
WEAPON  SYSTEMS  EQUIPMENT 

I 1A 


These  data  are  displayed  as  motion-time  histories,  shock  spectra, 
and  Fourier  magnitudes.  Two  examples  of  platform  motions  predicted  from 
measurements  will  be  discussed  in  this  section. 

Response  motions  were  determined  by  the  following: 

X^(t)  (•)  • T (oj)  X.  (w) 

where 

Xj^(t)  “ Accelerat ion- t ime  history  response  of  platform  at 
reference  point 

Xj^(<ii)  “ Acceleration-frequency  response  of  platform  at 
reference  point 

x.(u))  " Input  building  acceleration.  Annex  M or  local, 

Fourier- transformed  to  frequency  domain 

t((d)  ” Measured  overall  transfer  function  of  platform 

As  discussed  previously,  the  load  path  from  each  isolator  was  measured  from 
the  isolator/building  junction  to  a reference  point  on  the  platform.  These 
measurements  were  summed  for  all  isolators  of  a platform  to  give  an  overall 
(or  global)  transfer  function.  This  implies  that  the  building  input  motion 
is  uniform  and  in  phase  over  the  area  covered  by  the  platform.  Such  an 
assumption  was  required  because  the  building  areas  covered  by  each  platform 
specify  a uniform  input  into  each  isolator.  Alternatively,  each  transfer- 
function  path  and  the  building  input  motion  at  each  attachment  point  could 
have  been  Individually  computed  and  all  computations  summed  to  determine  the 
results.  The  dual  of  the  above  equation  is  the  convolution  integral  from 


which  the  identical  platform  response  may  also  be  computed  as: 


where 


t , T • Time 

X(t)  “ Acce lerat ion- t ime  history  response  of  platform  at 
reference  point 

x(ui)  “ Acce  lera  t ion- frequency  response  of  platform  at 
reference  point 

h(t)  “ Measured  impulse  function  of  platform 

(transformed  from  transfer  function  t(u)) 

X.(t)  “ Input  building  acceleration,  Annex  M or  local 

(Appendix  A) 

t “ Time  delay  function  for  the  convolution  operation 

In  effect,  this  convolution  corresponds  to  shifting  two  functions 
and  computing  the  integral  of  its  product  with  the  other  time  history  as  t 
ranges  from  0 to  T.  The  total  duration  of  platform  response  is  equal  to 
the  duration  of  the  impulse  function  plus  the  duration  of  input  time  history. 

Response  motion  predictions  of  the  control  room  platform  of  the  PAR 
Power  Plant  is  displayed  in  Figures  69  through  71.  Acceleration-time  histories 
for  the  frequency  band  35  Hi  to  500  Hi  and  0.5  Hi  to  500  Hi  are  included  in 
Figures  69a  and  69b.  In  the  displacement-time  history  (Fig.  69c),  it  is  noted 
that  the  peak-to-peak  displacement  is  less  than  +1  in.  These  peak  values  must 
be  considered  conservative  or  upper  bound,  in  view  of  the  non  I inear i t ies  in 
the  low  frequencies  as  previously  discussed.  Rigid-body  design  limits  for 
platforms  is  +6  in.  The  shock  spectrum  and  Fourier  spectrum  for  this  platform 
response  is  provided  in  Figures  70  and  71. 

Acceleration-time  history  response  of  Platform  D represents  two 
methods  of  determination.  Figure  72  was  developed  by  a finite-element  computer 
model  and  Figure  73  was  obtained  at  Waterways  Experiment  Station.  Both  time 
histories  eii.tjloyed  Annex  M input  motion  to  the  platform.  Rigid-body  modes 
and  properties  of  vertical,  roll,  and  pitch  were  analytically  inserted  for 
the  low  frequencies  of  the  measured  data  to  provide  comparability  between 
the  figures,  most  particularly  for  shock  spectra.  The  shock  spectra  of 
Figure  7^  are  computer-model  and  measured-data  response  predictions. 
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Platforn  PAftPP'CR*.  Predicted  ftcc*lerftt  ion-t  tm*  history  response,  using  th«  computed 
gtobel  transfer  function  convolved  with  Annex  H input  (frequency  range  3S  to  SOO  Hx) 
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(b)  Platform  PARPP-CR:  Accal era t ion- t ime  history  for  Annex  H Inpjt 
(frequency  range  O-S  to  SOO  Hx) 


(c)  Platform  PARPP-CR:  01 splacemenc*t ime  history  for  local  Input 
(frequoncy  range  0.5  to  SOO  Hz) 


'9.  RESPONSE  MOTION  PREDICTION  OF  CONTROL  ROOM  IN  PAR  POWER  PLANT 
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FIGURE  70.  PLATFORM  PARPP-CR:  VERTICAL-RESPONSE  SHOCK  SPECTRUM  FOR  ANNEX  M INPUT 
(FREQUENCY  RANGE  0.5  TO  500  HZ) 


0 so  100  ISO  200  2S0  SOO  3S0  MO  <S0  SOO 

ntQUCNCT,  HZ 

FIGURE  71.  PLATFORM  PARPP-CR:  PREDICTED  ACCELERATION  FOURIER  RESPONE  .MAGNITUDE, 
USING  THE  GLOBAL  TRANSFER  FUNCTION  CONVOLVED  WITH  ANNEX  M INPUT 
(FREQUENCY  RANGE  35  TO  500  HZ) 


I 18 


FIBURE  72.  PAR  POWER  PLANT,  PLATFORM  D:  ACCELERATION-TIME  HISTORY 
RESPONSE  COMPUTED  BY  FINITE  ELEMENT  METHOD 


FIGURE  73.  PAR  POWER  PLANT,  PLATFORM  D:  ACCELERATION-TIME  HISTORY 
RESPONSE,  ANALYTIC  0 TO  10  HZ  AND  TRANSFER  FUNCTION 
MEASUREMENT  10  HZ  TO  150  HZ 


FIGURE  JU.  PAR  POWER  PLANT,  PLATFORM  D:  FINITE-ELEMENT  SHOCK  SPECTRUM 

RESPONSE  CALCULATION  COMPARED  WITH  SHOCK  SPECTRUM  FROM  TRANSFER 
FUNCTION  MEASUREMENTS 
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Differences  result  from  the  waveforms  of  the  two  methods  and  the  higher 
magnitudes  of  motion  prediction  for  the  finite  element  case. 

2.2.3  RESPONSE  PREDICTION  Of  A SOIL-COVERED  ARCH  TO  AIR  BLAST 

This  unpublished  information  was  kindly  provided  by  Robert  E.  Walker 
of  the  U.S.  Army  Waterways  Experiment  Station,  Vicksburg,  Mississippi.  Several 
1/3-scale  model  hanqarettes  (Fig.  50)  were  subjected  to  air-blast  pressures  in 
EVENT  mixed  COMPANY  (500-ton  TNT  explosion).  Air-blast  pressure  records  taken 
from  one  of  the  hangarettes  are  given  in  Figures  75a  to  75d. 

Mobility  measurements  (V/F)  were  made  from  the  inside  surface  of 
the  reinforced-concrete  arch  to  soil  surface  near  the  locations  of  the 
air-blast  pressure  gages  on  a full-size  structure.  Using  these  measurements 
in  a reciprocity  sense  (assumed  linearity),  response  velocity  predictions 
were  made  from  the  mobility  measurements  taken  on  the  full-scale  structure. 

The  air-blast  data  was  scaled  to  prototype.  For  comparison  to  velocity  records 
of  the  arch  in  the  MIXED  COMPANY  Tests,  the  full-scale  predictions  were  scaled 
to  the  models.  The  scaled  response  velocity  at  the  mid-arch  position  is  given 
in  Figure  76.  To  provide  for  comparison  between  MIXED  COMPANY  recorded 
velocity  and  the  predictions  of  Figure  76  the  predicted  velocity-time  history 
had  to  be  lowpass  filtered  to  match  the  frequency  characteristics  of  the 
velocity  gages  used  in  the  test.  Test  data  and  predictions  are  shown  in 
Figures  77  and  78. 

Walker  used  the  response  prediction  procedure  set  forth  in  Equation  8 
of  Section  2.2.1,  except  for  the  local  engulfment  function.  This  omission  may 
account  for  the  late-time  oscillation  occurring  at  125  msec  in  Figures  76 
and  78.  The  prediction  using  mobility  measurements  in  addition  to  the  scaling 
employed  shows  surprising  good  agreement  with  actual  data.  It  can  also  be 
concluded  that  soil  nonlinearity  for  this  test  and  at  these  pressure  levels 
exerts  only  a moderate  influence.  Late-time  displacement  of  the  test  record 
(275  msec)  cannot  be  accounted  for  at  this  time,  since  it  is  not  known  whether 
the  gage  record  has  been  corrected  for  offset  and  drift.  Additionally,  the 
mobility  measurements  were  not  low  enough  in  frequency  pick-up  ground  roll. 

(Text  continued  on  p.  125-) 
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TIME,  MSEC 

FIGURE  77.  VELOCITY  GAGE  RECORD  AT  MID-ARCH  OF  HANGARETTE  FROM  MIXED  COMPANY  EVENT 
(500-TON  TNT  TEST) 
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Tint  IN  NUC 

FIGURE  78.  VELOCITY  MOTION  PREDICTION  FROM  MOBILITY  (V/F)  MEASUREMENTS 
ON  HANGARETTE  (Velocity  low-pass  filtered  to  match  velocity 
gage  frequency  response  of  Fig.  77)  (Courtesy  Waterways 
Experiment  Station) 
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2.3  SIMULATION  OF  STRUCTURAL  MOTIONS 


2.3.1  METHOD 

Transient  shock  tests  on  equipment  and  systems  to  simulate  the 
motions  induced  by  a nuciear  event  or  an  earthquake  are  largely  limited  to 

•k 

single-axis  test  machines.  Further  limitations  exist  in  the  size  and  weight 
of  equipment  which  can  be  tested.  Simulating  multiaxis  loading  on  large 
equipment  with  many  degrees  of  freedom  represents  a difficult  problem,  because 
it  is  impractical  to  generate  continuously  varying  forces  of  sufficient  mag- 
nitude. This  problem  is  further  aggravated  where  in-place  or  field  tests 
are  required.  On  the  other  hand,  short-duration  forces  of  large  magnitudes 
over  a wide  frequency  range  can  be  generated  by  mechanical  pulse  generators 
(Ref.  6).  Since  a discrete  number  of  pulses  superficially  presents  an 
appearance  quite  different  from  a continuous  excitation  signal,  it  becomes 
necessary  to  select  the  pulses  in  such  a way  that  the  resulting  vibration  of 
the  structure  matches  as  closely  as  possible  the  response  (e.g.,  displacement, 
velocity,  or  acceleration)  produced  by  the  continuous  force,  as  determined 
by  an  appropriate  error  criterion.  This  approach  is  shown  in  Figure  79. 

It  is  important  to  note  that  the  method  of  Figure  79  requires  that 
the  criterion  response  to  the  continuous  input  be  known,  which  would  generally 
not  be  true  in  practice.  To  accomplish  this  objective,  the  approach  proposed 
here  assumes  that:  (1)  a mathematical  model  of  the  system  under  study  is 
known,  and  (2)  the  inputs  of  interest  (e.g.,  earthquake  or  nuclear  blast)  are 
given.  Under  these  conditions  the  "criterion  response"  can  be  calculated 
and  used  to  obtain  the  pulse  train  for  the  simulated  test. 

In  general,  the  response-time  history  of  a test  article  under 
simulated  test  should  show  a reasonable  approximation  to  the  expected 
environmental  phenomena  for  meaningful  hardness/vulnerability  evaluation. 


A few  biaxial  and  triaxial  test  machines  are  available. 


2.3.2  PULSE-SIMULATION  ENVIRONMENTAL  TESTS 

The  need  for  a mechanical  Torce-pulse  generator  arose  in  connection 
with  studies  on  how  to  effectively  simulate,  by  physical  test,  shock  tran- 
sients on  massive  equipment  located  in  protective  structures.  The  force 
generators  could  not,  for  physical  reasons,  be  interpyosed  between  the  equip- 
ment and  the  building.  The  platform  environments  represent  the  building 
motions  generated  by  nuclear  attack  as  transmitted  into  the  equipment  through 
the  isolation  system. 

By  drawing  a metal  bar  or  mandrel  through  a cutting  tool  (or  vice 
versa)  with  suitable  motive  power  (air  pressure,  hydraulic  pressure,  explo- 
sive force,  electrical,  mechanical),  a series  or  a set  of  force-time  his- 
tories may  be  generated  (Refs.  6,  39)-  Reaction  at  the  attachment  points  of 
the  device  transmits  a force  output  to  the  structure  under  test.  Figure  80 
illustrates  the  device  using  motive  power  supplied  by  the  stored  energy  in 
a pneumatic  cylinder.  Amplitude,  duration,  and  shape  of  the  pulse  time 
histories  are  controlled  by  the  relative  velocity  between  the  cutting  tool 
and  the  metal  projection  on  the  mandrel,  and  by  the  shape  of  the  metal  pro- 
jection on  the  mandrel.  A photograph  of  the  pulse  generator  Installed  for 
in-place  testing  is  shown  in  Figure  81. 

Large  forces  may  be  generated  from  this  device.  Force  required  to 
cut  metal  is  largely  independent  of  rate  (velocity),  and  is  a function  of 
the  volume  of  chips  cut  (depth,  width,  and  length  of  cut)  dnd  a function  of 
the  specific  energy  of  cutting  a material.  A load  cell  or  strain  gage  may 
be  incorporated  in  series  with  the  device  to  provide  a force-time  history 
readout  as  the  device  is  operated. 

Another  form  of  pulse  generator  is  sketched  in  Figure  82  and 
pictured  in  Figures  83  and  81*.  This  form  is  a variation  of  the  drop  shock 
test  machine  but  differs  by  generating  multiple  pulses.  Expected  capacity 
of  the  design  of  Figure  82  is  70,000-lb  force;  and  with  design  modifications 
of  the  nubbins  and  cutter,  150,000-lb  force  output  is  expected.  Currently 
the  machine  is  used  to  generate  very  low  frequency  response  of  large  structures 

(Text  continued  on  p.  133.) 
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FIGURE  8k.  MEDIUM  FORCE  PULSE  GENERATOR 
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where  the  force  output  of  vibration  machines  is  very  low.  This  drop  shock 
generator  is  so  new  that  performance  information  is  limited  to  the  oscillo- 
graph data  of  Figure  85.  Figure  85  records  data  taken  on  the  Electrical 
Distribution  Center  shown  in  Figures  kS  and  46. 

Four  pulse  generators  of  Figure  8l  were  used  for  shock- i sol ated 
platform  tests  (Fig.  39)  and  were  attached  near  the  four  corners  of  each 
platform.  The  schematic  of  the  pulse  generator  and  power-actuation  system 
Is  given  in  Figure  86.  Power  to  the  pulse  generators  was  provided  by 
hydraulic  cylinders  having  volumetric  compensators  in  series.  All  four  units 
were  connected  in  series  by  metal  piping.  A pneumatic-hydraulic  accumulator 
provided  stored  energy  to  drive  the  generators  at  line  pressures  to  1000  ps i . 
Pu 1 se- in i t iat i on  timing  was  accounted  for  by  pre-positioning  each  cutter/ 
mandrel  and  performing  calibration  runs  of  the  entire  pulse  system  for  each 
platform  configuration. 

For  the  in-place  pulse-simulation  process,  a discrete  number  of 
pulses  appears  quite  different  from  a continuous  excitation  signal.  It  is 
necessary,  therefore,  to  select  the  pulses  in  such  a way  that  the  resulting 
vibration  of  the  platform  matches  as  closely  as  possible  the  response 
(i.e.,  displacement,  velocity,  or  acceleration)  produced  in  the  platform  by 
the  continuous  excitation  resulting  from  the  building  motion  (nuclear 
threat  input).  The  accuracy  of  simulation  is  determined  by  an  appropriate 
error  criterion  as  shown  in  Figure  79. 
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FIGURE  85.  OSCILLOGRAPH  RECORDS  AMPLITUDE-TIME  HISTORIES  OF  DROP  CARRIAGE  AND  EDC  STRUCTURE 
RESPONSE 
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The  basic  criterion  used  for  simulation  accuracy  is  the  integral 
squared  error  between  the  criterion  and  pulse-simulated  platform  response. 
The  error  function  is  evaluated  at  a sufficient  number  of  points  within  the 
mu  1 t i p 1 e-deqree-of - f reedom  platform  to  characterize  the  platform  as  com- 
pletely as  possible.  With  the  error  criterion  given,  the  pulse  occurrence 
time,  pulse  widths,  and  the  pulse  amplitudes  are  selected  by  a systematic- 
search  algorithm  such  that  the  error  is  minimized  (Ref.  40) . In  order  to 
use  the  simulation  method  in  conjunction  with  the  optimization  procedure 
discussed  above,  the  following  steps  were  performed; 

a.  The  impulse  functions  for  each  platform  were  determined  for 
each  pulse  location.  This  was  accomplished  by  converting 
measured  transfer  impedance  functions  in  the  frequency 
domain,  to  transfer  impedance  impulse  functions  in  the 
time  domain.  Typical  functions  for  each  platform  are  shown 
in  Figures  87  and  68. 

b.  Using  the  optimization  algorithm,  the  criterion  platform 

' response  was  converged  upon  by  pulse  trains  convolved  with 

the  above  impulse  function.  Typical  computed  pulse  trains 
are  shown  in  Figures  89  and  90. 

Using  the  pulse  profiles  specified  above,  piston  velocities  of  the 
hydraulic  system  (Fig.  86)  were  established  and  mandrels  of  the  pulses 
were  machined.  Measured  pulse  trains  from  in-place  tests  are  also  shown  in 
Figures  89  and  90-  Typical  acceleration-time  histories  obtained  are  dis- 
played in  Figures  91  and  92.  These  figures  are  a three-way  comparison  of 
(1)  criterion  (predicted)  response  to  nuclear  threat,  (2)  pulse-simulated 
' response,  and  (3)  actual  (measured)  response  of  platforms  to  in-place 

pulse  tests.  The  same  three-way  comparison  for  each  platform  is  also  shown 
in  shock  spectrum  format  in  Figure  93. 


(a)  Magnitude,  | x (u) /F  (ui)  ] 


(b)  Phase,  4>{u)) 


(c)  Impulse,  h [X(t)/F(t)l 


FIGURE  88.  CONTROL  ROOM  PLATFORM:  TYPICAL  TRANSFER  INERTANCE 
FUNCTION,  FROM  PULSER  ATTACHMENT  POINT  TO  REFERENCE 
ACCELERATION  POINT,  NO.  1,  ON  PLATFORM 
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SECTION  3 


IMPEDANCE  USAGE  AND  APPLICATIONS 

The  previous  section  set  forth  practical  engineering  projects  of 
impedance,  and  of  motion  predictions  and  environmental  simulation  based  on 
impedance.  A notable  feature  of  some  of  these  projects  was  their  dismayingly 
large  size  and  massiveness.  Numerous  data  measurements  were  necessary  for 
such  projects,  with  subsequent  complex  chains  of  multiplications,  divisions, 
and  summations  that  required  compatibility  in  both  the  time  and  the  frequency 
doma ins. 


However,  the  demonstrated  practicality  of  impedance  and  mobility 
provides  means  by  which  certain  types  of  problems  may  become  tractable.  For 
example,  the  results  of  previous  tests  could  be  processed  by  impedance/ 
mobility  for  scaling  the  data  from  one  type  of  ship  for  application  to  another 
type,  or  for  equipment  changes  and  replacements.  Other  applications  would 
include  quasi-l inear  systems,  nonlinear  systems,  parameter  identification, 
multivariate  regression  analysis,  and  pulsed  environmental  simulation.  These 
areas  are  further  discussed  in  this  section. 

Throughout  this  report,  the  term  impedance  has  been  used  in  a broad 
sense  because  it  is  widely  accepted  and  understood.  In  practice,  however, 
mobility  and  its  derivative  inertance  are  applied  almost  exclusively,  not 
only  for  ease  of  computation  and  for  visual  interpretation  but  also  because 
they  are  independent  of  the  number  and  locations  of  other  measurement  points 
(Ref.  Al). 

3. I LINEAR  SYSTEMS 

Section  2 covered  the  use  of  impedance  and  transfer  functions  for 
response-motion  prediction  based  on  linear  or  quasi-linear  systems.  Currently 
five  major  projects  on  this  subject  have  been  completed  or  are  in  the  process 
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I of  completion,  with  results  that  are  considered  quite  good.  As  a practical 

matter,  many  systems  exist  where  the  main  load-bearing  structure  (ship  or 
submarine  hulls)  remains  essentially  linear  over  a wide  range  of  attacks, 
yet  internal  equipment  may  exhibit  nonlinear  response  with  failure,  malfunc- 
tion, or  degradation.  The  nonlinear  interaction  of  equipment  on  the  linearly 
responding  main  load-bearing  structure  needs  to  be  assessed  on  a case-by-case 
basis.  Previous  projects  indicate  that  these  nonlinear  effects  are  not 
severe  enough  to  adversely  affect  impedance  measurements.  Response  predic- 
tions made  at  equipment  locations  at  criterion  or  postulated  attack  levels 
provide  input  functions  for  individual  labora tory-env i ronmen ta 1 tests  on 
the  equipment. 

Response  motions  (acceleration-time  histories)  determined  from 
impedance  measurements  may  be  made  quite  accurately  for  a large  class  of 
linear  and  quasi-l inear  systems.  The  accuracy  depends  on  tne  formulation 
of  appropriate  input  loads.  Scale-model  tests  for  air-blast  and  under- 
water shock  have  been  found  to  be  very  effective  for  formulating  the  input 
to  represent  the  transient-load  engulfment  of  the  target  facets.  Imprn"e- 
ment  in  the  accuracy  of  impedance  data  using  time-series  digital  data 
analysis,  noise  minimization  procedures,  and  new  measurement  devices  largely 
accounts  for  the  practicality  of  applications  in  predicting  response  motions. 
Each  succeeding  project  exhibits  continuing  improvements  in  reduction  of 
costs,  ease  of  testing,  and  accuracy.  In  the  near  future,  many  analytic 
models  will  be  checked  against  impedance  measurements  for  accuracy. 

3.2  NONLINEAR  SYSTEMS 

New  methods  of  approach  can  be  employed  that  provide  a means  of 
improved  predictions  for  nonlinear  response  of  structures.  By  proper  design 
of  tests,  data  can  be  extracted  in  a specific  manner  to  yield  nonlinear 
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describing  functions.  These  functions,  which  are  discussed  below,  quantita- 
tively specify  inelastic  and  damage  phenomena  and  permit  nonlinear  response 
predict  ions. 

In  the  parametric  ident i f icat i on  approach,  the  mathematical  struc- 
ture of  the  model  is  postulated  but  its  parameters  are  not.  Most  of  the 
work  done  so  far  in  st rue tura I - system  identification  has  been  considered  from 
the  parameter  estimation  approach.  The  identification  task  in  this  parametric 
model  approach  eventually  reduces  to  a search  in  parameter  space  where  system 
parameters  are  iterated  repeatedly  until  values  are  obtained  that  meet  a 
specified  error  criterion.  Among  the  techniques  used  in  this  approach  are 
gradient  or  random  search  methods. 

In  the  describing  function  identification  approach,  no  a priori 
assumptions  are  made  regarding  the  structural  configuration  of  the  mathemati- 
cal model  of  the  system  to  be  identified.  The  identification  problem  becomes 
a search  in  "functional  space."  The  basic  mathematical  statement  supporting 
this  process  is  presented  below. 

Let  u(t)  be  a random  noise  input  signal,  and  g.  be  a function  of 
n variables.  The  response  y(t)  of  a nonlinear  system  can  be  expressed  as 

n 

y (t)  “ Z)  G. [g. , u(t) 1 

I » I 

where  G.[g.,  u(t))  is  a functional  of  the  i degree,  which  is  orthogonal 
to  all  functions  G.,  j < i: 


G. [g. , u(t)]  X Gj(q..  u(t)]  ” 0 for  i < j 


I'*? 


The  first  few  terms  of  the  functionals  are: 


G,(g,.  u(t)J  ^y*  g,(r,)  u(t  - t^)  dr^ 


G^lgj,  u(t)  ] 


’ fj  ' ''2^ 


- 92^'’2'^2^  ^^2 

ijii 

Gj[g^,  u(t)]  “ y^yy9j(t  J .T^  ,T  j)  u(t  - t^)  u(t  - t^)  U(t  - T^)  dTjdt2dT. 

- 3iyJ'gj(t,  .T2.T2)  u(t  - T|)  dtjdT^ 


where  k 

t 


average  power  of  random  noise. 


Since  the  parametric  identification  requires  a priori  assumptions 
regarding  the  characteristics  of  the  system,  large  errors  can  be  introduced 
in  modeling  complex  structures  if  the  order  of  the  mathematical  model  does 
not  agree  with  the  order  of  the  actual  structure.  This  problem  becomes 
particularly  acute  for  the  response  in  the  nonlinear  range.  However,  for 
the  describing  function  identification  approach,  the  order  of  the  system 
need  not  be  assumed,  but  is  allowed  to  develop  to  optimally  satisfy  the 
measured  characteristics. 


For  the  general  solutions  of  nonlinear  systems,  the  Volterra  Series 
may  be  used  to  represent  an  explicit  input/output  relation.  One  of  the  main 
problems  in  the  application  of  the  theory  is  the  explicit  determination  of 
the  kernels  appearing  in  the  Volterra  Series.  A solution  improvement  set 
forth  by  Wiener  was  to  expand  these  kernels  in  terms  of  a set  of  orthogonal 
functions,  such  as  Laguerre  functions.  Practical  solution  was  achieved  by 
cross  correlation. 
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Studies  of  nonlinear  effects  (Ref.  1)  disclosed  non  1 inear i t ies 
that  could  be  very  sensitively  detected  by  rapid  sine  sweeps  (CHIRP  tests) 
without  requiring  excessively  large  forces  that  could  drive  the  test  structure 
into  the  damage  region.  From  this  information,  an  alternative  procedure  has 
evolved  in  the  form  of  state  variables,  by  which  Chebychev  polynomials  are 
used  as  functionals  relating  input  to  output. 

3.3  PARAMETRIC  IDENTIFICATION 

By  extending  parameter  identification  to  nonlinear  systems,  a series 
of  impedance  measurements  at  different  force  levels  is  obtained.  This 
approach  provides  a group  of  nonparamet r i c impedance  curves  from  which  parame- 
ters are  identified.  From  the  parameters  obtained,  functional  relations  are 
set  forth  in  the  model,  as  for  nonlinear  springs  and  damping,  to  provide  an 
approximate  model  of  the  nonlinear  transmission  path.  The  parametric  approach 
can  also  be  used  in  conjunction  with  the  functional  approach  to  determine 
system  parameters  once  the  functionals  (describing  functions)  are  determined. 
This  combined  approach  presages  an  improved  means  of  extracting  physical 
phenomena  for  application  to  nonlinear  modeling  of  structures. 

Random  search  algorithms  for  parameter  optimization  have  been  widely 
applied  and  documented.  They  have  the  advantage  of  (1)  leading  to  global 
solutions  of  nonlinear  systems,  (2)  guaranteeing  convergence,  and  (3)  easing 
computer  implementation.  On  the  negative  side,  random  search  algorithms  may 
converge  very  slowly,  particularly  in  criterion  surfaces  of  high  dimensionality. 
Several  procedures  have  been  tried  in  the  past  to  circumvent  slow  convergence. 

The  algorithm  used  in  several  impedance  and  pulse  simulation  projects 
is  another  approach  to  the  determination  of  the  optimal  step  size  (Ref.  40). 
Rather  than  a fixed-length  step,  steps  used  are  random  in  both  length  and 
direction.  Hence,  the  adaptation  described  below  is  based  on  the  selection 
of  the  optimal  variance  of  the  step-size  distribution  as  the  search  progresses. 
Large  variances  are  desirable  in  the  early,  exploratory  portions  of  a search. 
However,  in  the  vicinity  of  a local  optimum,  a smaller  value  of  the  standard 
deviation,  o,  will  decrease  the  probability  of  overshoot. 
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The  algorithm  for  the  adaptive  random  search  consists  of  alternating 
sequences  of  a global  random  search  with  a fixed  value  for  the  step  size 
variance  followed  by  searches  for  the  locally  optimal  o.  This  adaptive 
random  search  method  was  used,  as  shown  in  the  following  example,  to  optimize 
17  parameters  of  a nonlinear  soil  model. 


In  specific  applications,  data  are  given  in  the  form  shown  in 
Table  where  c and  o are  the  vertical  soil  strain  and  stress,  respectively, 
obtained  from  a uniaxial  test,  and  where 


a.  ,0 . 
' J 


First  stress  invariant  (lo.) 

Second  reduced  stress  (stress  deviation) 

\ 2 
invariant  “ -r  Z [a.  - a.) 

b I j 

Principal  stress  components 


The  stress  invariants  depend  on  e,  o,  and  the  loading  history,  as 
well  as  numerous  other  soil  parameters  that  are  to  be  identified.  Among  the 
present  soil  models  that  are  being  widely  used  in  ground-shock  computations 
are  elastic,  perfectly  plastic,  variable  moduli,  and  an  elastic-plastic  model 
with  a movable  cap . 


Once  a candidate  mathematical  model  is  selected,  the  problem  is 
then  to  determine  the  "best"  set  of  parameters  that  will  result  in  a loading/ 
unloading  curve  and  stress  path  that  will  simultaneously  match  both  sets  of 
experimental  data  such  as  given  in  Table  4.  Thus,  the  identification  problem 
is  reduced  to  a search  in  parameter  space  that  will  yield  the  extremum 
(minimum)  value  of  a suitable  penalty  function. 

The  cap  model  used  to  represent  soils  employs  a yield  surface  that 
combines  both  strain  hardening  and  ideal  plasticity.  The  ideally  plastic 
modified  Orucker-Prager  criterion  represents  the  ultimate  shear  strength  of 
the  material  and  is  associated  with  fracture  or  sustained  plastic  flow  in 
laboratory  experiments.  The  form  of  the  yield  criterion  is 

f,  - (J,,  Jp  - 0 
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TABLE  A.  REFERENCE  DATA 


Strain 

1 ncrement 

1 

' Strain 

Stress 

'1 

-1 

.000- 

03 

r 

.000-03 

-1 

.669-03 

-3. 

,000-03 

5. 

77A-OA 

-6 

.900- 

02 

! -7 

.000-02 

-2 

.097-02 

-3. 

890-02 

6. 

,928-03 

3 

.000- 

03 

-6 

.700-02 

-5 

.000-03 

-9. 

,000-03 

1 , 

.732-03 

-5 

.000- 

03 

-7 

.200-02 

-1 

.828-02 

-3. 

,290-02 

6, 

.371-03 

-2 

.800- 

02 

-1 

.000-01 

-3 

.8AI-O2 

-7. 

,A90-02 

1 , 

.155-02 

-5 

.000- 

02 

-1 

.500-01 

-7 

.552-02 

-1 . 

,500-01 

2, 

. 19A-02 

-2 

.500- 

02 

-1 

.750-01 

-1 

.080-01 

-2. 

,2A9-01 

2, 

,887-02 

-2 

.300- 

02 

-1 

.980-01 

-1 

. AOO-01 

-3. 

,000-01 

3. 

, A6A-02 

-1 

.200- 

02 

-2 

. 100-01 

-1 

.680-01 

-3. 

,7A9-01 

3. 

,811-02 

-1 

.500- 

02 

-2 

.250-01 

-1 

.960-01 

-A. 

500-01 

A. 

,0Al-02 

-5 

.000- 

03 

-2 

. 300-01 

-2 

.220-01 

-5. 

,2A9-01 

A. 

.157-02 

-5 

.000- 

03 

-2 

.350-01 

-2 

.A9O-OI 

-6. 

,000-01 

A, 

.328-02 

-5 

.000- 

03 

-2 

.AOO-01 

-2 

.760-01 

-6. 

,7A9-01 

A. 

AAA-02 

-5 

.000- 

03 

-2 

.A50-01 

-3 

.010-01 

-7. 

.500-01 

A. 

.503-02 

-1 

.000- 

03 

-2 

.A6O-OI 

-3 

.270-01 

-8, 

,2A9-01 

A. 

.563-02 

-1 

.000- 

03 

-2 

.A7O-OI 

-3 

.596-01 

-9. 

.180-01 

A, 

.619-02 

2 

.000- 

OA 

-2 

.A68-01 

-3 

. 2A0-01 

-8, 

.2A9-OI 

A. 

.272-02 

2 

.000- 

OA 

-2 

.A66-01 

-2 

.950-01 

-7. 

.500-01 

3. 

.982-02 

3 

.000- 

OA 

-2 

.A63-01 

-2 

.670-01 

-6. 

.7A9-OI 

3. 

.635-02 

3 

.000- 

OA 

-2 

.A60-01 

-2 

.389-01 

-6. 

,000-01 

3 

.AOA-02 

1 

.000- 

03 

-2 

.A5O-OI 

-2 

. 110-01 

-5, 

,2A9-01 

3. 

.118.02 

1 

.000- 

03 

-2 

.AAO-01 

-1 

.830-01 

-A, 

,500-01 

2, 

,387-02 

1 

.000- 

03 

-2 

.A30-01 

-1 

.5A0-01 

-3. 

,7A9-01 

2, 

,5A0-02 

.000- 

03 

-2 

.A20-01 

-1 

.260-01 

-3. 

,000-01 

2, 

.250-02 

.000- 

03 

-2 

.AlO-Ol 

-9 

.683-02 

-2, 

.2A9-01 

1 . 

.903-02 

5 

.000- 

OA 

-2 

.AO5-OI 

-6 

.910-02 

-1  , 

,500-01 

1 . 

,672-02 

5 

.000- 

OA 

-2 

.AOO-01 

-A 

.083-02 

-7. 

,A90-02 

1 , 

,386-02 
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As  the  plastic  volumetric  strain  increases  or  decreases,  the  strain-hardening 
cap  expands  or  contracts,  respectively.  The  cap  is  denoted  by 


(J,,  J'.  c ) 


where  c is  plastic  strain. 

It  is  assumed  that  elastic  behavior  is  governed  by 


do.  . = Kde,  I 6.  . + 2G  (de.  . - -5-  dc.  .6.  .] 

I J kk  I J I J 3 kk  I J 


where 


do.  . 
I J 

de . 


I J 


dc 


kk 


6.  . 

I J 

K 

G 


Stress  increment  tensor 

Strain  increment  tensor 

Linear  dilatation  (=  de.|^  + *^^22  ^ '^'"33^ 

Kronecker  delta  (“1  if  i » j,  » 0,  if  i j) 

Bulk  modulus,  which  is  a function  of  pressure 

Constant  shear  modulus 


The  basic  cap  model  used  in  this  study  consists  of  the  following 
parts:  (1)  a variable  bulk  modulus,  (2)  a constant  shear  modulus,  (3)  a 

fracture  surface,  and  (A)  a cap.  The  cap  model  is  defined  by  the  following 
equations  and  an  associated  flow  rule: 


Bulk  and  Shear  Modulus 


where 


Blj’- 


Be'^*^’  BeV’l 
®2^  - 


G 

B^  through  B^ 


Constant 

Empirical  coefficients 
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Fracture  Surface  ( f , ) 


n/T- 


Y,  - V^e 


Cap  (f  ) 


(L  - J,)^  Y^  - (L  - x)^ 


- ^)  - 


L > Y^\Y2e  - YJ  - X = 0 


,P  P P P 

'l  ' * ^2  ^ ^3 


where 


First  invariant  of  plastic  strain 

'S'^  P 


Z-  e 
i-1 


1 • 


Distance  to  the  center  of  the  elliptical 
cap 


X = Intersection  of  elliptical  cap  with 

J ^ axis 

Y^  through  Y^  • Empirical  coefficients 

In  the  application  of  this  model  to  soils,  the  cap  movement  is 
controlled  by  the  change  of  plastic  volumetric  strain;  hence,  strain  hardening 
can  be  reversed. 
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The  computer  implementation  of  this  model  results  in  a code  that, 
when  supplied  with  a strain  time  history,  will  generate  the  corresponding 
state  of  stress,  and  subsequently  the  stress  invariants,  at  every  point  in 
the  stress  path. 


In  view  of  the  complexity  of  the  aforementioned  model, 
implementation  has  required  a substantial  effort  and  resulted  in 
computer  code  that  performs  a significant  number  of  computations 
incremental  change  in  the  strain  level  in  order  to  determine  the 
state  of  stress. 


its  computer 
a large 
for  every 
correspond ing 


Given  n discrete  data  points  with  coordinates 


[e. , Oj (e.  , ?)]  and 


J,  (Ci,  ?). 


of  the  type  shown  in  Table  A where 


c . 

I 

? 

o,(e,,  ?) 
J,  (e.,  ?) 
J-  (Ci,  ?) 


Value  of  independent  variable  (strain),  i-l,  2,...,  n 

Parameter  vector  whose  k components  correspond  to 
the  parameters  to  be  identified 

Digitized  values  of  stress,  i“l  , 2,...,  n 

Digitized  values  of  stress  invariant 

Digitized  values  of  stress  invariant 


The  penalty  function  J used  in  the  present  example  is: 


where 


S,  - [o(E.)-o(e.))^ 

I I I 


S2  = i:  [J,(e,)  - J,{e.)1^ 


Wj , W^,  = Weighting  functions 

n = Number  of  data  points 


and  where  c,  o,  represent  measurement  data  and  o,  are  the 

corresponding  quantities  obtained  by  specifying  the  sequence  of  strain 
increments  Ac.  = e.  , - e.,  i - 1,2,...,  n - 1,  and  then  determining  the 
state  of  stress  in  accordance  with  the  predictions  of  a suitable  mathematical 
mode  I . 

The  random  search  algorithm  was  used  to  fit  soil  data  by  performing 
the  following  steps: 

a.  The  discretized  values  of  the  stress/strain  and  the  stress 
path  of  Table  k were  used  as  data. 

b.  A judicious  choice  was  made  for  the  soil  cap  mode)  fitting 
parameter  values,  which  were  subsequently  used  as  initial 
parameter  values  to  start  the  random  search  algorithm.  These 
initial  values  are  listed  in  Table  5. 
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c.  Using  the  initial  values  given  in  Table  5 in  conjunction  with 
an  appropriate  cap  model  computer  code,  Figure  9^*  shows  a 
comparison  of  this  "initial  solution"  and  the  experimental 
data.  Note  that  this  "initial  solution"  was  obtained  by 
merely  substituting  the  parameter  values  given  in  Table  5 into 
the  material  package  subprograms.  The  optimization  section  of 
the  parameter  identification  code  was  not  exercised  in  this 
step. 


d. 


Using  equal  error  weighting  functions  Wj  = w^  = w^  = I in 
the  criterion  function,  the  deviation  error  between  the 
"initial  solution"  and  the  given  data  was  0.272. 


e.  Again  using  the  "initial  values"  listed  in  Table  5 as  a starting 
estimate  to  initiate  the  parameter  identification  code  and  then 
operating  the  code  for  a few  hundred  iterations,  an  optimum 
set  of  parameters  was  found.  These  parameters  are  listed  in 
Column  3 of  Table  5-  The  optimized  solution  is  shown  in 
Figure  95  and  its  corresponding  deviation  error  is  0.070 
{'  17k  of  the  starting  error). 

Note  that  the  reduction  of  the  deviation  error  was  achieved  by 
optimizing  the  values  of  17  parameters  listed  in  Table  5-  The 
percentage  change  in  these  parameters  ranged  from  -12.8% 
to  +13%. 


3 ■ k SCALING  AND  CORRELATION 

3.k.l  MULTIVARIATE  REGRESSION  ANALYSIS 

Multivariate  regression  analysis  may  be  viewed  as  the  problem  of 
determining  a linear  relationship  between  two  or  more  variables.  It  may  also 
be  viewed  as  a linear  curve  fitting  procedure  in  which-- 

Y - + 6,Xj  + BjXj  + . . . + B^X^ 
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FIGURE  S‘t.  COMPARISON  OF  EXPERIMENTAL  DATA  AND  INITIAL  SOLUTION 


where 


Y ” Dependent  variables 

8 » Regression  coefficients  to  be  determined 

n 

X » Independent  variables 

n 

In  the  regression  computation,  values  of  6 are  assumed;  and  with 
the  independent  variables  of  X^,  Y is  computed.  This  Y is  compared  to 
the  known  Y value  for  the  platform  and  g values  are  adjusted  through 
iterative  procedures  to  within  some  error  limit.  The  independent  variables, 
X^,  are  the  particular  characteristics  of  any  structural  system. 

There  are  two  difficulties  that  arise  from  the  use  of  multivariate 
regression  analysis; 

• Proper  choice  of  the  model 

• Variability  of  the  results 

The  results  may  be  meaningless  if  an  improper  model  of  the  situation 
is  employed,  or  at  least  be  less  meaningful  than  some  other  model.  Stepwise 
regression  procedures  can  be  used,  and  a general  engineering  awareness  of  the 
d v.mensi  onal  i ty  of  the  problem  tends  to  overcome  this  first  difficulty. 

The  second  problem  Is  a function  primarily  of  the  number  of  struc- 
tural systems  used  in  the  computations:  the  larger  the  number  or  more 
uniformly  distributed  the  data  base,  the  better  the  results.  Procedures  for 
the  analysis  of  variability  of  errors  can  be  applied  to  judge  the  reliability 
of  the  results.  In  particular,  confidence  bounds  may  be  calculated  and 
placed  on  either  side  of  the  results.  These  confidence  bounds  are  usually 
given  at  a certain  level.  For  example,  one  frequent  way  of  calculating  them 
is  such  that  chances  are  the  true  results  will  be  inside  the  confidence 
bounds  95%  of  the  time. 
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Comparison  of  systems  is  made  both  by  similarities  and  by  differences. 
In  formulating  trial  models,  the  independent  parameters  reflect  the  similari- 
ties, and  the  specific  values  assigned  thereto  provide  the  differences.  In 
shock  or  Fourier  response- spec t ra-est i mat i on  procedures  for  a ship,  multi- 
variate regression  analysis  is  applied  once  for  each  frequency  being  examined. 
For  example,  the  X and  Y variables  could  be  defined  as  follows: 

Y “ Fourier  pressure  spectra  for  frequency  f.  on  ship's  hull 

Xj  = Depth-charge  size 

X^  * Depth-charge  depth  below  surface 

Xj  * Depth-charge  distance  from  snip 

X^  = Length  of  ship  (water  line) 

X^  = Beam  of  ship 

Xg  = Length/beam  ratio 

Xj  = Hull  t h i ckness 

Xg  = Displacement  vol/weight 

The  units  of  these  parameters  are  immaterial;  the  multivariate 

regression  analysis  procedure  automatically  accounts  for  the  units.  It  should 

be  noted  that  the  X variables  may  be  composite.  That  Is,  they  might  be 

2 3 2 A 

equal  to  functions  such  as  t + w , w f , exp(-af) , etc.,  where  4 is 
length,  w is  weight,  f is  frequency,  and  a is  a constant. 

Grouping  of  parameters  as  in  the  foiegoing  permits  formulation  of 
reasonable  representative  system  models  or  functional  relationships,  rather 
than  brute  force  polynomial  fits  to  data. 

For  example,  shock- i so  I a ted  platforms  protecting  equipment  comprise 
62  different  configurations  totaling  107  platforms  distributed  throughout  a 
hardened  facility  (Ref.  9,  30).  Prediction  of  response-motion  for  each 
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platform  and  for  each  environmental  location  within  the  system  was  required 
for  hardness/survivability  assessments  of  equipment.  The  a priori  selection 
of  platforms  as  a representative  sample  for  transfer-function  measurements 
permitted  interpolation  methods  to  be  used  to  estimate  with  acceptable  confi- 
dence the  response  motions  of  all  platform;.  A simple  environmental  scaling 
rule  then  allowed  these  estimates  of  platform  motions  to  be  adjusted  to  the 
threat  environments  for  a particular  building  and  floor  location. 

A multivariate  regression  method  was  used  to  predict  responses  for 
all  types  of  platforms,  measured  and  unmeasured.  The  platforms  can  be 
classified  in  terms  of  their  physical  characteristics  (i.e.,  length,  width, 
density,  etc.)  and  local  input  environments.  Transfer-function  measurements 
and  subsequent  response  analysis  of  the  10  selected  platforms  resulted  in 
10  shock- response  spectra  for  prediction  analysis. 

The  polynomial  approach  was  generated  by  representing  the  shock- 
isolated  platforms  in  the  form  of  plates.  Since  the  dynamic  response  of  a 
plate  (hence  its  shock  spectrum)  is  directly  related  to  its  frequency  spectrum, 
it  is  reasonable  to  choose  the  following  characteristic  parameters  as  inde- 
pendent parameters  for  the  regression  analysis:  t,  w,  o,  wt/iso.  On  the 
basis  of  several  trial  cases,  the  following  nonlinear  multivariate  regression 
was  evolved; 


x(ss)  “ 6q  + 8,1  + 82^^  + 8,4^  + 8,^0  + 6r0^ 


+ 8^P^  + 8^(4/w)  + 60(wt/lso) 


where 

x(ss) 

8 

n 

4 

P 

w 

wt/i so 


Acceleration  shock  spectra 

Coefficients  of  regression  for  each  frequency  point 
Length 

Weight/area  (calculated  density) 

Width 

Total  weight  per  isolator 
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The  p.irticulcir  data  bases  used  for  this  analysis  were  the  dependent 
variable  values  of  the  shock- spec t rum  responses  of  the  10  measured  platforms 
over  the  frequency  range  of  35  to  500  Hz,  and  the  independent  variable  plat- 
form characteristics:  length,  width,  density,  length/width,  and  weight  per 
isolator.  Character i st ics  of  the  measured  platforms  are  provided  in  Table  6. 
To  ensure  optimum  predictive  results,  the  error  limit  for  the  response  data 
was  set  at  +1.’  about  the  mean  for  the  10  measured  platforms.  Also,  because 
of  the  optimum  distributive  sampling  of  the  platforms  for  testing,  the  upt sr 
and  lower  bounds  of  the  62  sets  of  platforms  did  not  exceed  the  upper  and 
lower  bounds  of  the  10  platforms  investigated  (see  Fig.  96). 

Figures  97  and  98  show  two  platform  responses  to  the  referent  over- 
all environment  input  as  computed  from  impedance  and  transfer  function  data. 
Overlaid  on  each  plot  is  the  predicted  shock  spectrum  for  that  same  platform, 
derived  by  multivariate  regression  analysis. 

The  predicted  shock  spectra  show  a good  correlation  with  the  actual 
shock  spectra  if  all  10  platforms  are  viewed  in  a statistical  relationship. 
Individually,  7 of  the  10  platforms  are  almost  perfectly  predicted  to  within 
i1/1<  standard  deviations. 

The  shock-spectrum  responses  of  the  remaining  unmeasured  52  types 
of  platforms  were  predicted  by  using  their  specific  characteristics  (length, 
width,  length/width,  density,  wt/iso)  with  the  coefficients  generated  by 
multivariate  regression  analysis.  These  values  are  the  platform  responses  to 
the  overall  input.  Typical  examples  of  responses  predicted  for  unmeasured 
platforms  by  the  regression  method  are  shown  in  Figures  99  and  100. 

To  develop  the  platform  shock  spectrum  for  its  particular  location 
(local  shock  spectrum),  a scaling  technique  was  developed.  This  scaling 
method  was  used  to  develop  the  response  of  the  62  types  of  platforms  o the 
local  input  environments.  The  procedure  calls  for  dividing  the  shock  spectrum 
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FIGURE  96.  DISTRIBUTION  OF  SHOCK  SPECTRA  FOR  10  PLATFORMS  MEASURED  IN  PLACE, 
Response  computations  determined  from  overall  Annex  M input. 
Environmental  levels  for  equipment  tested  on  vibration  shakers 
are  also  shown  (Annexes  E,  H,  L,  »nd  K) . 


FIGURE  97.  PLATFORM  PARPP-O:  COMPARISON  OF  SHOCK  SPECTRA  DETERMINED  FROM  IMPEDANCE 
MEASUREMENTS  AND  FROM  REGRESSION-PREDICTION  METHOD 


m 

FIGURE  98.  PLATFORM  PARPP-CR:  COMPARISON  OF  SHOCK  SPECTRA  DETERMINED  FROM  IMPEDANCE 
MEASUREMENTS  AND  FROM  REGRESSION  PREDICTION  METHOD 
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FIGURE  99.  PLATFORM  MSCB-250A:  VERTICAL  SHOCK-RESPONSE  SPECTRA 
(REGRESSION  PREDICTION  METHOD) 
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FIGURE  100.  platform  MSRPP- I : VERTICAL  SHOCK-RESPONSE  SPECTRA 
(REGRESSION  PREDICTION  METHOD) 
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values  of  the  local  environments  by  the  shock  spectrum  values  of  the  overall 
(reference)  environment.  A relationship  can  be  drawn  as  follows  to  form  a 
pseudotransfer  function: 

"(ss)  "(ss)  local  input  ss 

'local  ref.  ref.  input  ss 

where 

= Platform  response  shock  spectrum  to  local  input 
* environment 

y^^*^  = Platform  response  shock  spectrum  to  overall 

(reference)  input  environment 

local  input  ss  = Input  shock  spectra  specified  for  a particular 
bu i I d i ng  and  f 1 oor 

ref.  input  ss  = Overall  (reference)  shock  spectrum 

Figures  99  and  100  also  show  platform  response  to  local  shock 
spectrum  input  environments  determined  from  the  scaling  techniques.  Scaling 
of  shock  spectra  as  above  must  be  performed  with  some  caution.  Where  large 
changes  in  waveform  and  time  durations  exist,  errors  in  response  estimates 
can  result.  However,  for  estimating  the  high-frequency  response  for  any 
isolated  platform,  this  is  a very  useful  and  low-cost  tool. 

3.^4. 2 IMPEDANCE  SCALING 

Systems  with  similar  dimensions  may  be  scaled  for  response  by  matching 
impedance  functions,  even  though  extensive  not-to-scale  elements  may  be 
present  (Ref.  32,  42). 

One  method  of  optimal  scaling  overlays  each  matched  pair  of  impedance 
» functions  from  naval  vessels  and  compares  them  for  the  best  curve  fit  (minimum 

error).  The  resulting  magnitude-scale  factor  and  frequency-scale  factor  are 
then  used  to  scale,  for  example,  such  data  as  the  response  due  to  a depth 
charge  on  one  submarine  to  predict  the  response  of  another  type  or  class  of 
- submarine. 
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Another  method  conforms  more  to  conventional  network  theory  in 
accordance  with  the  following  representation; 


Xj  (t)  (“)Xj  (ju>) 


f , (ju) 


where 


X 


1 


(t) 


X,  (jw) 

.i/j  ( ju) 
Fj  (jw) 


Acceleration  time  history  at  an  internal  location  in 
vessel 


Complex  acceleration-frequency  response  of  x(t) 


Generalized  inertance  functions 
measured  on  Vessel  No.  1 


( impedance) 


Generalized  input  force  on  “essel  No.  I 


For  the  second  vessel,  generalized  inertance  functions  (impedance) 

are  measured  and  formulated  such  that  response  prediction  on  the  second  vessel 

may  be  scaled  from  the  first  vessel  as: 


x^{t)  (»)  X2{ju)) 


,y/-(ju) 

.rf7n=T  h'J"' 


The  predicted  response  for  Vessel  No.  2 is  based  on  equivalent 
input  loading  for  Vessel  No.  1.  Depending  on  the  geometric  differences  of 
the  vessel,  scaling  may  also  be  required  for  input  loads.  Scaling  of  input 
loads  may  also  be  used  to  determine  equivalent  response  motion  (shock  or 
Fourier  spectra)  of  the  second  vessel.  The  inverse  problem,  i.e.,  defining 
the  external  loads  from  the  above  formulation,  is  quite  difficult  and  will 
not  be  looked  into  here. 

The  generalized  inertance  function  may  be  viewed  as  a network 
equivalent  system,  a matrix  array,  or  a sample.  The  form  of  representation  can 
be  developed  for  the  application  required.  In  References  32  and  A2,  single 
representative  samples  for  each  of  the  four  locations  scaled  were  used  with 
quite  good  results  even  though  their  factors  differed  significantly  (see 
Sec.  2.1.6. I , Table  2). 
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3.5  PULSE  APPLICATION 


Pulse  train  excitation  of  structures  was  covered  in  Section  1.1 
and  in  more  detail  in  Section  2.3.2.  In  addition  to  mechanical  pulse  genera- 
tors that  require  fixed  reaction  points  or  inertial  reaction  masses,  fluid 
and  chemical  pulse  generators  are  under  application  development  and  should  be 
available  in  the  next  12  to  15  months.  These  new  pulse  generators  depend 
on  the  mass-flow  reaction,  and  as  a consequence  permit  considerable  flexibility 
in  applications,  including  underwater  use.  Applications  will  determine  the 
pulse  generator  type  best  suited  for  use,  such  as: 

a.  Peak  force  required 

b.  Pulse  durations  (max  and  min) 

c.  Number  of  pulses  required 

d.  Turn-on/turn-of f timing 

e.  Test  bandwidth  frequency 

f.  Number  of  pulse  units  required 

g.  Number  of  test  axes  for  simultaneous  testing 

force  generation  may  be  configured  for  environmental  response 
simulation,  systems  (impedance)  measurements,  modal  surveys,  and  diagnostic 
testing.  Diagnostic  testing  is  performed  in  terms  of  structure  response,  of 
equivalent  sine,  chirp,  and  random.  The  latter  two  types  of  excitation  are 
particularly  suited  for  studies  of  nonlinear  phenomena. 

In-place  testing  of  large  equipment,  subsystems,  and  systems  can 
be  accomplished  by  pulse-simulation  methods.  Transmission  of  shock  loads 
from  the  hull  of  a vessel  through  the  structure  to  the  multiple-mounting 
points  of  a ship's  equipment  alters  the  input  to  complicated  acceleration 
time-history  vectors.  Section  k.]  covers  projects  for  predicting  these  time- 
history  motions  for  each  mounting  point  and  for  each  orthogonal  direction 
as  well  as  for  selected  points  on  the  equipment  itself.  Pulse  generators 
would  be  used  to  match  these  motion  predictions  in  a vector  sense  on  the 
equipment.  Physical  duplication  of  expected  or  actual  motions  of  equipment 
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provides  the  most  reliable  means  for  determining  failure,  malfunction, 
degradation,  and  damage  as  a function  of  attack  threat.  Experience  has  also 
shown  that  once  fragility  data  on  equipment  is  obtained,  hardness  upgrade 
may  be  significantly  improved.  Typical  nitty-gritty  details  of  a fragility 
test  and  upgrade  hardening  on  electronic  communication  equipment  is  given 
in  Tables  7 and  8 from  Reference  3-  This  transient  shock-fragility  program 
objectively  demonstrated  that  the  equipment  would  meet,  and  exceed  with  high 
probability,  the  weapon-system  mission  requirements  for  nuclear-weapon 
effects  of  shock  and  vibration. 

Fragility  of  systems/components  (Ref.  2)  is  categorized  in  terms 
of  failure,  malfunction,  and  damage/degradation.  These  terms  demonstrate  that 
fragility  is  more  comprehensive  than  expected  after  a cursory  throught:  its 
operational  connotation  extends  the  meaning  well  beyond  structural  breakage. 
However,  even  the  fragility  categories  are  subject  to  different  interpretation 
in  the  range  of  technical  fields  involved.  Therefore,  the  categories  are 
defined  as  fol lows: 

• Failure  is  defined  as  an  irreversible,  env i ronment- induced , 
inoperative  condition  or  operation  outside  of  tolerances. 

J rrevereible  refers  to  the  system/ component  remaining  inopera- 
tive or  out  of  tolerance  after  the  environment  is  removed. 

• Malfunction  also  represents  an  envi ronment- i nduced  inoperative 
condition  or  out-of-tolerance  operation.  However,  the  process 
is  reversible,  i.e.,  the  system/component  returns  to  satis- 
factory operation  upon  removal  of  the  environment. 

• Damage/degradation  has  multiple  meanings.  It  may  be  considered 
as  a mild  form  of  failure  (i.e.,  irreversible  but  borderline 
operation).  It  also  represents  permanent  degradation  in  per- 
formance, reduction  in  hardness,  or  limitation  of  survivability. 
Damage  also  applies  to  system/ component  attributes  that  are 
unrelated  to  performance. 
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TABLE  7-  qualification  (1004)  LEVEL  TEST  SUMMARY  (Ref.  3) 
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TABLE  8.  FRAGILITY  (200%)  LEVEL  TEST  SUMMARY 
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In  matching  the  response  of  a structure  to  a specific  input-dynamic 
environment,  it  is  recognized  that  this  approach  is  not  as  fruitful  for 
nonlinear  systems  as  it  is  for  linear  systems.  It  is  planned  to  use  the 
pulse  technique  to  compare  calculated  {from  nonlinear-mathematical  models) 
and  measured  responses  due  to  applied  diagnost ic-pui se  trains.  The  results 
would  be  approximate  with  the  resulting  error  as  an  index  to  the  degree  o*' 
nonlinearity  encountered.  Thus,  a base  is  provided  (if  the  approxima t i on 
is  poor)  for  improving  the  analysis  and  validating  or  identifying  the 
nonlinear  system  characteristics. 
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SECTION  k 

APPLICATIONS  TO  ONGOING  NAVAL  NEEDS 

Intensive  and  extensive  applications  of  impedance/mobility  can  be 
made  to  ongoing  naval  needs.  The  benefits  that  justify  such  efforts  are 
Outlined  briefly  in  the  nine  projects  suggested  in  this  section. 

Although  impedance  has  a solid  analytic  foundation  (see  p.  6),  the 
examples  shown  in  Section  2 and  the  projects  described  in  this  section  are  all 
based  upon  physical  measurements.  Impedance  data  provide  useful  knowledge 
of  systems  in  as-built  configurations.  In  contrast,  analytical  models  must 
use  engineering  material-properties  data,  estimates  of  joints  and  workmanship, 
and  modeling  simplification.  The  need  to  use  abstract  or  simplified  data 
rather  than  actual  measurements  has  given  rise  in  recent  years  to  probabilistic 
solutions  with  associated  confidence  levels.  The  large  uncertainties  that 
result  may  lead  to  pessimism  for  hardness  and  survivability.  An  interesting 
check  on  the  validity  of  probabilistic  solutions  could  be  the  collection  and 
statistical  analysis  of  matched  impedance  measurements  from  ships  of  the 
same  class  and  type. 

Finite  element  models  and  impedance  methods  can  be  complementary. 

For  the  same  bandwidth,  impedance  measurements  may  be  used  to  verify  and 
improve  the  finite  element  models.  For  high  frequencies,  impedance  measure- 
ments are  more  effective;  whereas  for  highly  nonlinear  systems,  finite  element 
methods  are  preferable.  The  two  methods  are  not  mutually  exclusive  but  can 
be  melded  together  in  the  frequency  bands  where  each  is  most  effective. 

We  feel  that  important  information  can  be  gleaned  from  the  following 
projects.  They  are  suggested  for  review,  extension,  and  elaboration  by  the 
naval  community. 


k. \ IMPEDANCE  APPLICATIONS 

k.\.\  PROJECT  1;  PREDICTION  OF  BASE  (ATTACHMENT-POINT)  MOTIONS  OF  THE  SHIP'S 
STRUCTURE  WHEN  LOADED  WITH  EQUIPMENT  OR  FOR  CHANGE  OF  EQUIPMENT 

a.  Purpose : To  determine  Input  motion  to  equipment  for  devising 
shock  and  vibration  tests  on  this  equipment  for  vulnerability 
and  hardness  requirements. 

b.  Theory : The  shock  environment  at  equipment  locations  in  a 
naval  vessel  is  specified  from  tests  and  procedures.  Dynamic 
interaction  at  the  mounting  locations  is  modified  by  interaction 
of  the  ship's  structure  and  equipment.  If  new  equipment  is 

to  be  substituted  and  dynamically  differs  from  the  original 
equipment,  then  the  originally  specified  motion  is  altered. 

c.  Procedure : Inertance  measurements  are  made  on  the  ship's 
structure  and  on  the  equipment.  Ratios  of  these  functions 

are  multiplied  by  the  free  (unloaded)  velocity  (or  acceleration) 
motion  of  the  ship  at  the  mounting  locations  to  predict  the 
motion  with  equipment  attached. 


v'(t)(-)  V'(u) 


V(u.) 


('i-D 


where 

v'(t)  • Velocity  time  history  at  ship's  structure/ 

equipment  mounting 

V (li))  ■ Velocity-frequency  spectrum  of  v'(t) 

V(u)  > Veioci ty-f requency  spectrum  of  ship  structure 
only 
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'*^22^'^^  = Driving-point  inertance  of  ship's  structure 

' Driving-point  inertance  of  equipment 

d.  Demonstration  and  Verification:  Demonstration  and  verifica- 
tion of  this  project  may  be  performed  on  the  U.S.  Navy  Medium- 
Weight  Shock  Machine  and  on  the  U.S.  Navy  Floating  Shock 
Platform  (FSP)  or  the  Large  Floating  Shock  Platform  (LFSP). 
Tests  would  be  made  with  and  without  equipment  mounted  on 
both  the  shock  machine  and  the  FSP/LFSP,  and  velocity  records 
would  be  taken.  Subsequent  inertance  measurements  would  be 
used  in  Equation  A-1  to  predict  motions.  These  predictions, 
when  compared  to  actual  test  records,  would  determine  the 
accuracy  and  suitability  of  the  method. 

e . Applications  to  the  Fleet: 

• Standardized  inertance  measurement  procedures  and  equip- 
ment could  be  established  to  determine  threat  motions 

to  equipment. 

• The  ensemble  of  threat  motions  on  equipment  could  be  used 
to  establish  survivability,  vulnerability,  and  fragility 
by  tests  in  naval  laboratories  or  by  naval  contractors. 

• Hardening  and  upgrading  of  marginal  equipment  may  be 
i mp I emen ted. 

• Shock  and  vibration  transmitted  from  ship's  structure  to 
equipment  can  be  attenuated  by  the  use  of  isolation  and 
damping  devices. 

• Shock-sensitive  equipment  may  be  relocated  to  locations 
more  environmentally  benign  in  a ship. 
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PROJECT  2:  SCALING  HULL  PRESSURE  AND  RESPONSE  OF  NAVAL  VESSELS  TO 
EXPLOSIVE  LOADINGS 

a.  Purpose.  To  scale  the  response  motions  and  hull  pressure 
loadings  generated  by  explosive  loads  on  one  type  of  naval 
vessel  to  similar  or  different  types  of  vessels,  when  not-to- 
scale  effects  exist. 

b.  Theory.  In  addition  to  presently  available  scaling  procedures, 
inertance  seal ing  and  multivariate  regression  analysis  may  be 
used  on  an  empirical  basis. 

c.  Procedure:  Procedures  for  multivariate  regression  analysis 
and  inertance  scaling  are  given  in  Section  3- 

d.  Demonstration  and  Verification:  An  array  of  inertance 
measurements  taken  on  two  different  ships  are  matched  for 
approximately  equivalent  paths  (e.g.,  control  room  to  similar 
hull  position).  Ratios  of  these  paired  measurements  are  used 
to  scale  the  response  of  one  vessel  to  the  other.  The 
predicted  responses  are  compared  for  accuracy  to  actual  test 
records . 

e.  Applications  to  Fleet:  Empirical  procedures  will  be  developed 
and  verified  by  impedance  methods  and/or  multivariate 
regression  analysis  to  scale  explosive  tests  on  naval  vessels 
to  various  classes  and  types. 

Multivariate  regression  analysis  requires  a reasonable  data  base 
and  has  the  characteristics  of  correlation.  Inertance  scaling  is  functionally 
related  and  requires  only  measurements  of  the  two  systems  that  are  to  be 
re lated. 
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Multivariate  regression  analysis  requires  a data  base  sufficient  to 
make  reasonably  accurate  predictions-  It  may  be  more  appropriate  to  examine 
hull  response  or  pressure  loading  as  the  dependent  variable,  since  inertance 
measurements  and  scaling  can  be  used  to  determine  responses  elsewhere  within 
the  ship.  Parameters  required  would  include  charge  weight,  depth  of  charge 
(below  surface),  target  distance,  depth  of  water,  target  orientation  (including 
surface  or  below  surface),  hull  parameters,  and  vessel  velocity.  Iterative 
models  of  the  parameters  would  be  generated,  and  would  be  expected  to  be  a 
combination  of  linear  and  nonlinear  terms. 

Of  particular  importance  is  the  scalability  of  the  shock  motions  of 
equipment  tested  on  the  FSP  or  LFSP  to  their  locations  on  naval  vessels.  In 
another  sense,  for  the  same  or  dynamically  similar  equipment  motion,  what  is 
the  explosive  threat  required  to  induce  this  motion  on  another  ship? 

if. 1.3  PROJECT  3:  PREDICTION  BY  IMPEDANCE  TECHNIQUES  OF  NAVAL-VESSEL 
RESPONSE  TO  EXPLOSIVE  CHARGES 

a.  Purpose : To  determine  acceleration  time  histories  and 
associated  acceleration  Fourier  frequency  spectra  and  shock 
spectra  at  critical  equipment  locations  in  U.S.  naval  ships 
and  submarines.  These  responses  can  be  determined  for  a range 
of  threats  and  for  various  orientations  of  the  vessel  with 
respect  to  the  threat  source.  Both  single  attack  and  multiple 
attacks  from  several  directions  can  be  predicted.  Respjonse 
predictions  would  be  based  upon  quasl-l inear  hull  response. 

b.  Theory : Response  motions  at  selected  interior  points  of  a 
naval  vessel  may  be  predicted  by  complex  multiplication  of 
transfer  impedance  functions  and  transient  hydraulic-pressure 
loadings  over  the  surface  of  the  hull.  Transient-pressure 
loadings  on  the  hull  may  be  formulated  analytically  from 
scale  model  tests  and  from  low-level  depth  ch.  'ges  on  full- 
size  naval  vessels. 
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P rocedure : Transfer  impedance  measurements  are  taken  from 
critical  equipment  location  to  the  ship's  hull,  as  illustrated 
in  Figure  101.  Spacing  of  the  measurements  is  determined  by 
the  gradient  observed  between  contiguous  measurements  on  the 
hull  and  the  engulfing  velocity  of  the  shock  wave  from  the 
weapon  (see  Fig.  58).  The  impedance  functions  are  arranged 
with  respect  to  shock-wave  engulfment  of  the  vessel.  It  is 
to  be  noted  that  the  surface  area  of  the  hull  of  a Polaris  sub- 
marine has  about  the  same  area  as  the  roof  of  the  PARE  building, 
which  was  measured  for  impedance/motion  prediction  (see 
Sec.  2. 2. 1.2).  Hull  pressure-t ime  histories  for  a range  of 
attack  levels  and  directions  would  be  established  at  input 
threat  criteria  and  obtained  from  scale  model  tests,  calcula- 
tions, and  full-scale  tests. 

d.  Demonstration  and  Verification:  Demonstration  and  verification 
of  the  procedures  may  be  applied  to  any  naval  vessel  that  was 
used  for  explosive  tests  and  for  which  test  records  exist  or 
for  which  new  tests  can  be  made.  As  a practical  convenience, 
the  Floating  Shock  Platform  (FSP),  the  Large  Floating  Shock 
Platform  (LFSP) , and  the  Submersible  Submarine  Test  Vehicle 
(SSTV)  can  be  used.  The  resulting  acceleration  time-history 
records  from  a test  would  be  compared  to  predicted  motions 
based  on  pressure  loadings  and  impedance  functions  to  assess 
the  degree  of  accuracy  obtained. 

e.  Application  to  Fleet: 

• Impedance  measurements  on  naval  vessels  can  be  made  with 
minimal  disruptions  of  operations  and  no  damage/risk 
even  for  measurements  at  high-input-force  levels. 

• Response  predictions  can  be  made  prior  to  explosive  tests 
on  surface  vessels  and  submarines  to  assess  the  safety 

of  the  test. 
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• Worse-case  single-  and  multiple-attack  scenarios  can  be 
developed,  and  response  levels  can  be  compared  to 
laboratory-determined  damage  levels  of  equipment  for 
hardness/survi vabi 1 i ty  assessments. 

• Response  predictions  may  be  scaled  to  other  naval  vessels 
as  discussed  in  Project  2. 

• Improvements  can  result  in  test  criteria  and  environmental 
specifications  for  shipboard  equipment. 

Projects  1 through  3 are  predicated  on  linear  and  quasi-l inear  sys- 
tems. "Quas i - 1 i near"  is  taken  in  the  sense  that  the  ship's  hull  remains  linear, 
and  other  elements  in  the  system  may  respond  partially  into  the  nonlinear  range, 
as  in  damping  and  stiffness.  For  cases  where  non  1 inear  1 1 ies  are  greater  than 
known,  response  motions  predicted  will  be  conservative  (larger)  with  some  fre- 
quency shifts  expected.  This  situation  will  predominate  in  early  times. 

4.1.4  PROJECT  4:  NONLINEAR  SYSTEMS 

a.  Purpose : To  extract  functionals  relating  input/output  for  sub- 
system's components  or  equipment  exhibiting  large  non  1 i near i t i es . 

b.  Theory:  Several  methods  exist  for  the  development  of  func- 
tionals, such  as  Wiener  kernels  and  state  variables.  These 
can  be  combined  with  multilevel  testing  to  generate  nonlinear 
parameter  models  from  measured  data. 

c.  Procedure : General  procedures  are  as  covered  in  Section  3- 

d.  Demonstration  and  Verification:  A project  has  not  yet  been  com- 
pletely defined  for  this  subject.  However,  demonstration  and 
verification  can  be  applied  to  the  medium-weight  shock  machine 
when  driven  to  its  stops  and  to  shipboard  piping  systems. 

e.  Applications  to  Fleet: 

• Nonlinear  shock- i solat ion  systems  for  shipboard  equipment 

• Shipboard  piping  systems 
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4.2  SHOCK-SIMULATION  PROJECTS 


4.2.1  PROJECT  5;  pulse-simulation  TESTS  ON  U.S.  NAVY  MEDIUM-WEIGHT  SHOCK 
MACHINE 

a.  Purpose : To  demonstrate  pu I se-s imu I a t i on  methods  and  to 
increase  the  testing  flexibility  of  the  medium-weight  shock 
machi ne. 

b.  Procedure : Response  of  the  test  structure  and  anvil  is 
calculated  for  a pressure  time  history  as  shown  in  Figure  102. 

A pulse  train  is  computed  to  generate  an  equivalent  response 
from  which  pulse  tests  are  conducted. 

c.  Ver i f i cat i on : Accuracy  of  the  test  is  determined  by  comparing 
the  predicted  motions  to  test  records. 

4.2.2  PROJECT  6:  PULSE-SIMULATION  TESTS  ON  MAJOR  SUBSYSTEMS 

a.  Purpose;  To  qualify  a major  subsystem  to  criteria  shock  loads, 
to  determine  hardness/fragility  levels,  to  identify  failure 
modes,  and  to  upgrade  hardness  levels  where  required. 

b.  Procedure : On  selection  of  the  major  subsystem  for  testing-- 
for  example,  the  missile-launch  tube  shown  in  Figure  103-- 
response  motions  determined  from  the  previous  impedance  studies 
of  Project  4 are  used  to  generate  pulse-equivalent  motions. 
Shock  tests  are  performed  simultaneously  in  the  three  ortho- 
gonal axes.  Testing  may  be  accomplished  as  installed  aboard 
ship  if  space  permits,  or  at  a Navy  yard,  naval  laboratory,  or 
contractor  facility. 

c.  Ver i f icat ion : The  tests  are  validated  by  comparison  of  test 
motion  to  prediction  and  to  environmental  tests  on  the  vessel. 
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Pulse-driven  anvil  and 
test  structure 


FIGURE  102.  PULSE  DRIVEN  MED lUM-WE I GHT  SHOCK  MACHINE  SIMULATING  UNDERWATER 
SHOCK  MOTION  ON  TEST  STRUCTURE 
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k.Z.}  PROJECT  7:  PULSE-S I MULAT I ON  TESTS  ON  PIPING  SYSTEMS 


a.  Purpose : To  qualify  piping  systems  to  criterion  shock  loads, 
to  determine  hardness/fragility  levels,  to  identify  failure 
modes,  and  to  upgrade  hardness  levels  where  required. 

b.  Procedure : As  in  previous  tests,  pulse  methods  will  be  used 

to  duplicate  responses.  However,  for  piping  systems,  responses 
must  be  generated  from  the  quasi-l inear  to  strongly  nonlinear 
regions.  Both  nonlinear  models  and  nonlinear  functionals 
(from  test  measurements)  will  be  used.  Shock  testing  would 
be  performed  in  a naval  laboratory  or  shipyard  on  a piping 
system,  as  illustrated  in  Figure  lOA,  that  duplicates  a shipboard 
installation  for  safety  and  test  efficiency. 

c.  Val idat ion:  Validation  for  the  levels  of  shock  severity 
involved  could  best  be  made  by  comparing  the  pulsed-test 
piping-system  responses  to  responses  obtained  on  the  Large 
Floating  Shock  Platform  (lFSP). 

A. 2. A PROJECT  8:  PULSE-S I MULAT 1 ON  TESTS  ON  SHIPBOARD  RADAR  ANTENNA 

a.  Purpose : To  qualify  shipboard  radar-antenna  system  to  criterion 
air-blast  loads,  to  determine  hardness/f rag i I i ty  levels,  to 
identify  failure  modes,  and  to  upgrade  hardness  levels  where 
requ i red. 

b.  Procedure:  Response  prediction  and  pulse  simulation  of 
antenna  systems  are  as  depicted  in  Figure  105  and  require 
both  impedance  measurements  and  air-blast  pressure  load- 
distribution  information.  Air-blast  load  data  are  best 
obtained  from  scale-model  tests  from  the  large  shock  tube  at 
Dahlgren,  Virginia,  and  supplemented  by  full-scale  tests 

in  such  high-explosive  events  as  DIAL  PACK  and  DICE  THROW. 
Experimental  air-blast  data  on  truss-type  structures  will 
provide  necessary  information  on  transient-drag  coefficients 
and  shading  effects. 
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ELEVATION  VIEW 


FIGURE  104.  TESTS  ON  SHIPBOARD  PIPING  SYSTEM 
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FIGURE  105.  IN-PLACE  TESTS  OF  MAJOR  SUBSYSTEMS 
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k.l.s  PROJECT 

a. 


b. 


Va 1 i dat i on : Validation  of  pulse-simulated  responses  would  be 
limited  to  scale-model  tests  in  shock  tubes  and  to  full-scale 
tests  in  high-explosive  events. 

9:  PULSE-SIMULATION  TESTS  ON  A WEAPON'S  SYSTEM 

Purpose:  To  simulate  underwater  shock  response  of  a naval 
vessel . 

Procedure : Response-motion  predictions  as  performed  in 
Project  3 ere  used  to  specify  the  pulse  generators  and  their 
spatial  distribution  over  the  ship's  hull  (Fig.  106a).  For 
manned  or  tethered-operat iona I vessels,  pulse  testing  would  be 
at  a level  to  preclude  damage  consistent  with  safety.  High- 
level  tests  would  be  accomplished  on  the  Submersible  Submarine 
Test  Vehicle  (SSTV) , Floating  Shock  Platform  (FSP),  and  Large 
Floating  Shock  Platform  (LFSP)  Pressure- load  distribution 
over  the  hull,  as  shown  in  Figure  106b,  must  be  determined 
from  analysis,  scale-model  tests,  and  full-scale  depth  charges. 

Val idation:  Validation  of  pulse-simulated  responses  is  made 
by  comparison  to  equivalent  depth-charge  tests  on  the  weapon's 
system  or  test  vehicle. 
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PULSE  SIMULATION  TO  GENERATE  STRUCTURAL  RESPONSE  MOTIONS  EQUIVALENT  TO 
MOTIONS  INDUCED  BY  NUCLEAR  OR  DEPTH  CHARGE  UNDERWATER  SHOCK 


DETIRMINE  AND  SPECIFY  INTERFACE  PRESSURE  LOADS 


FIGURE  106.  WEAPONS  SYSTEMS  TEST 


SECTION  5 


RECOMMENDATIONS 

This  report  has  summarized  current  activities  in  the  measurement  of 
structural  dynamic  functions,  practical  methods  of  response  prediction  to 
environmental  threats,  and  a new  means  of  environmental  testing.  The  scaling 
of  systems,  particularly  by  impedance/mobility  techniques,  also  offers 
opportunities  for  refining  information  from  ongoing  and  past  projects  for 
hardness  assessment  and  upgrade. 

It  is  shown  in  this  report  that  mechanical  impedance/mobility  tech- 
niques have  come  of  age.  Impedance  and  mobility,  long  a tool  of  the  experi- 
menter, have  evolved  to  practical  methods  under  engineering  conditions.  The 
ability  to  accurately  predict  and  simulate  accelerat ion- time  histories  at 
the  mounting  locations  of  weapons  system  equipments  and  components  provides 
additional  knowledge  for  protecting  and  hardening  these  elements  to  meet 
mission  requirements  of  naval  ships  and  submarines. 

It  is  therefore  recommended  that  the  responsible  authorities  in  the 
Navy  Codes  consider  projects  for  applying  the  metliods  presented  and  for 
elaboration  of  them  for  specific  naval  projects. 

The  projects  outlined  in  the  previous  section  started  with 
demonstration/proof  examples  and  proceeded  to  hardness/survivability  levels 
for  naval  subsystems  and  systems.  The  following  recommended  projects  are 
listed  by  priority,  starting  with  the  less  complex  tests  to  ensure  proper 
application  of  the  methods,  procedures,  and  techniques  involved.  The  next 
group  provides  for  measurements  and  tests  where  a potentially  high  yield  of 
information  may  be  obtained  for  a minimal  expenditure  of  funds.  The  final, 
more  complex  projects  then  benefit  from  evolving  refinements  in  test  and 
measurements  and  the  learning  curve  of  personnel  and  of  management. 
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(1)  U.S.  NAVY  MEDIUM-WEIGHT  SHOCK  MACHINE  (PROJECTS  1 AND  5) 

• Determine  motion  of  shock  machine  as  altered  by  addition 
of  equipment  mounted  thereon. 

• Pulse-test  the  shock  machine  with  mounted  equipment,  to  simulate 
response  motions  induced  by  an  underwater  depth  charge. 

(2)  FLOATING  SHOCK  PLATFORM/LARGE  FLOATING  SHOCK  PLATFORM  (PROJECT  1) 

• Determine  motions  at  equipment  mounting  points  on  FSP/LFSP 
with  and  without  mounted  equipment. 

(3)  SCALING  RESPONSE  MOTIONS  (PROJECT  2) 

• Scale  response  motions  from  one  naval  vessel  to  another  and 
verify,  by  explosive  tests  on  second  vessel,  that  the  predicted 
motions  are  accurate  within  acceptable  error  bounds. 

• Scale  transient-shock  tests  made  of  equipment  on  the  Floating 
Shock  Platform  or  Large  Floating  Shock  Platform  to  shipboard 
installations  and  verify  with  test  records. 

• Scale  explosive  tests  for  Floating  Shock  Platform  or  Large 
Floating  Shock  Platform  to  yield  the  equivalent  motion  (energy 
spectra,  Fourier  spectra,  shock  spectra)  of  equipment  aboard 

a ship. 

{k)  RESPONSE  PREDICTIONS  TO  UNDERWATER  DEPTH  CHARGE (S)  (PROJECT  3) 

• Convolve  pressure  loads  from  depth  charge(s)  with 
systems  (impedance)  functions  of  naval  vessels  to  predict 
motions.  Candidate  vessels  are  submarine.  Submersible 
Submarine  Test  Vehicle  (SSTV),  Floating  Shock  Platform,  Large 
Floating  Shock  Platform,  or  other  vessels  to  be  designated. 
Predicted  motions  are  to  be  verified  by  tests. 
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NONLINEAR  RESPONSE  (PROJECTS  k AND  7) 

• Predict  motions  of  shipboard  piping  system  by  extraction  from 
data  of  nonlinear  functionals.  Pulse  simulation  of  nonlinear 
motion  of  shipboard  piping  system  would  duplicate  a shipbo  rd 
installation  and  be  tested  in  a laboratory. 


PULSE  SIMULATION  TESTS  OF  MAJOR  SUBSYSTEMS  (PROJECT  6) 

• Make  simultaneous  tests  in  three  orthogonal  axes  by 
pulse-simulated  response-motion  equivalents  on  major  sub- 
systems in  the  weight  range  of  20,000  lb  to  1(00,000  lb.  These 
tests  may  be  accomplished  in  place  as  installed  aboard  ship, 
in  naval  laboratories,  naval  shipyards,  or  at  contractors' 
plants.  Verification  of  simulation  will  be  by  comparison  to 
motion-response  records  from  shock-barge  tests. 

RADAR  ANTENNA  (PROJECT  8) 

• Perform  pulse-simulation  tests  on  deck-nounted  equipment  of 
surface  ships  by  simultaneous  tests  in  three  orthogonal 
axes.  These  tests  Induce  responses  caused  by  air-blast 
loading  and  may  be  performed  in  place  and  as  installed. 

PULSE  SIMULATION  OF  COMPLETE  WEAPONS  SYSTEMS  (PROJECT  9) 

• Test  a surface  ship  or  submarine  to  equivalent  underwater 
depth  charge(s)  by  attachment  of  pulse  generators  to  the 
hull.  Repeat  tests  for  different  levels  of  attack,  for 
different  directions,  and  for  multiple  attacks.  For 
environmental  protection,  this  test  method  would  kill  fewer 
fish  than  present  methods. 
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